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Based on a functional-integral formalism, a generalization of the self-energy-functional theory 
(SFT) is proposed which is applicable to systems of interacting electrons with disorder. Similar 
to the pure case without disorder, a variational principle is set up which gives the physical (dis- 
order) self-energy as a stationary point of the (averaged) grand potential. Although the resulting 
self-energy functional turns out to be more complicated, the formal structure of the theory can 
be retained since the unknown part of the functional is universal. This allows to construct non- 
perturbative and thermodynamically consistent approximations via searching for a stationary point 
on a restricted domain of the functional. The theory and the possible approximations are worked 
out for models with local interactions and local disorder. This results in a derivation of different 
mean-field approaches and various cluster extensions, including well-known concepts as the statisti- 
cal dynamical mean-field theory, the molecular coherent-potential approximation and the dynamical 
cluster approximation. Due to the common formal framework provided by the SFT, one achieves 
a general systematization of dynamical approaches, i.e. approaches based on the spectrum of one- 
particle excitations. New mean-field and new cluster schemes naturally appear in this framework 
and complement the existing ones. Their prospects for future applications are discussed. 
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I. INTRODUCTION 

The combined influence of electron-electron interaction 
and of disorder on material properties represents a cen- 
tral question of solid-state theory. In diluted magnetic 
semiconductors^ like Gai_ x Mn x As or Zni_j;Mn x Se, the 
magnetic properties, such as the Curie temperature, sen- 
sitively depend on the (random) distribution of the Mn 
ions as well as on the type and the strength of their ef- 
fective magnetic interaction which results from strong 
Coulomb interaction among the Mn 3d valence elec- 
trons. Several transition-metal oxides^ with partially 
filled metal 3d shells (e.g. manganites or cuprates) are 
antifcrromagnctic Mott or charge-transfer insulators and 
exhibit a rich phase diagram upon doping with charge 
carriers. The disorder potential introduced due to the 
substitution process considerably affects their magnetic, 
charge and orbital ordering. 

While these examples show the need for a comprehen- 
sive theory of interacting and disordered electron sys- 
tems, they also demonstrate the immense complexity one 
faces in any theoretical approach. Even strongly simpli- 
fied (Anderson-Hubbard-type) models with local inter- 
actions and local disorder only, arc highly non-trivial if 
studied in a regime which excludes a simple perturbative 
treatment .MiiLSiZ 

For three (and higher) dimensions, one may in first 
place focus on the local charge and spin dynamics of 
the electrons and, complementary to scaling theories^ 3 - 
disregard the long-wavelength modes which govern the 
immediate vicinity of a phase transition. In this con- 
text, mean- field approaches and cluster extensions are 
well justified. The mean-field concept is formally valid in 
the limit of high dimensions. Subsequent cluster exten- 
sions are suited to reincorporate short-range correlations 



which are neglected in the purely local mean-field ap- 
proach. Clearly, a mean-field treatment excludes impor- 
tant effects such as the destruction of long-range order 
due to thermal order-parameter fluctuations^ or Ander- 
son localization^ for example. Nevertheless, tractable 
mean-field theories can be valuable tools for an under- 
standing of interacting and disordered systems with dif- 
ferent competing orders and complex phase diagrams. 

A mean-field theory can be formulated on the level 
of Hamiltonians and electronic states. This yields sim- 
ple approaches such as the Hartree-Fock appoximation to 
treat the interaction part and the virtual-crystal approx- 
imation to treat the disorder part of the problem. 10 ! 11 
These are completely static theories which in addition to 
spatial also neglect temporal fluctuations. 

Temporal degrees of freedom can be taken into ac- 
count in a mean-field theory when this is based on the 
spectrum of excitations. Placing the one-particle Green's 
function in the center of interest, results in a mean-field 
theory which is distinguished by the fact that it yields 
the exact result in the limit of infinite spatial dimen- 
sions. With a proper scaling of the model parameters 
this limit preserves a highly non-trivial dynamiesi 12 ! 13 
This distinguished mean-field theory, for the interaction 
part of the problem, is the dynamical mean-field theory 
(DMFT) JdiiSiiSiii it gi V e S the exact (local) interaction 
self-energy of the prototypical Hubbard mode l 18 d 9 i 20 in 
the D = oo limit. For the disorder part, the coherent- 
potential approximation (CPA ) 21 ' 22 i 23 ' 24 gives the ex- 
act (local) disorder self-energy of the disorder Anderson 
model* in the D = oo limit. 

Phenomena depending on dimension are missed in a lo- 
cal mean-field approach but can be restored step by step 
using cluster expansions ! 11 ^ 5 Using a single-site mean- 
field theory as a starting point for a systematic expan- 
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sion is surely an inadequate approach to include long- 
wavelength modes and their effects. The main motiva- 
tion for the subsequent inclusion of spatial correlations in 
cluster theories is rather the expected rapid convergence 
of local observables. 

The purpose of the present paper is to contribute to a 
systematization of mean-held approaches and their clus- 
ter extensions in the combined case of interactions and 
disorder and to explore new approximation schemes. The 
strategy is to seek for a proper generalization of the self- 
energy-functional theory (SFT) developed recently . 26 > 27 
For the pure (disorder-free) case, it has been shown 
that different mean-held and cluster approaches are re- 
covered, and new approximations can be constructed 
in a systematic way which guarantees thermodynami- 
cal consistenc y 28 ' 29 Here, we describe a novel derivation 
of the SFT which is non-perturbative, i.e. a formulation 
which does not refer to formal sums of skeleton diagrams. 
This formulation is well suited for an extension of the 
theory to disordered (and interacting) systems. The gen- 
eralized SFT is worked out in detail. It is shown that it 
makes contact with (i.e. rederives) a number of previous 
approaches: 

(i) the DMFT+CPA put forward by Janis and 
Vollhardt^ - and by Dobrosavljevic and Kotliar , 31 ! 32 
which has recently been used to study metal-insulator 
transitions^ at non-integer filling and the effects of dis- 
order on magnetism ) 34 ' 35 ' 36 

(ii) the local distribution approach of Abou-Chacra, 
Anderson and Thouless^ which has recently been evalu- 
ated numerically by Alvermann and Fehske^ and which 
is the conceptual basis for 

(hi) the statistical DMFT proposed by Dobrosavljevic 
and Kotlia r 39 ' 40 with several recent applications,— e.g. 
to strongly coupled disordered electron-phonon systems, 

(iv) the molecular CPA 1 - 1 - and its combination with 
the cellular DMFT (C-DMFT) 4 ^ 3 . of Kotliar et al. and 
Lichtenstein and Katsnclson, 

(v) the dynamical cluster approximation (DCA) for 
disordered systems as introduced by Jarrcll, Krishna- 
murthy and Maier , 44 ' 45 

(vi) the disorder analog of a simplified DCA recently 
introduced by Minh-Tieni 6 - 

as well as with several variants of these approaches - 
such as the typical medium theory (TMT) 4 ^ involving 
the geometrical averaging of the local density of states 
which has been suggested by Dobrosavljevic, Pastor and 
Nikolic and applied in combination with DMFT by By- 
czuk, Hofstetter and Vollhardt.i 8 -^ 9 - 

The construction of a generalized SFT provides a uni- 
fied theoretical framework which is able to rederive and 
thereby to classify the above-mentioned approximations 
(i)-(iv) . This procedure automatically discloses the view 
on new approximations: Generalizations of the peri- 
odized cellular DMFT (PC-DMFT)^ and of the cluster- 
perturbation theory (CPT) 51 i 52 as well as the variational 
cluster approach (VCA) 28 i 53 i 54 are suggested for disor- 
dered (and interacting) systems. 



The main intention of the paper is to work out the 
formal concepts. The benchmarking and application of 
the different approaches requires a numerical implemen- 
tation which is beyond the present scope but intended 
for the future. 

The paper is organized as follows: The next Sec. |TT] 
introduces a number of basic quantities needed for the 
subsequent construction of the self-energy functionals. 
In Sec. IHII the self-energy functional for pure systems is 
derived non-pcrturbativcly within the functional-integral 
formalism. A brief general discussion of approximations 
follows in Sec. IIVI for the case of a fixed disorder config- 
uration. This provides the basis for the statistical SFT 
(statSFT) in Sec.|V]and for the statistical DMFT in par- 
ticular. The main ideas for the construction of the gen- 
eralized self-energy functional of configuration indepen- 
dent self-energies are provided in Sec. I VII while Sec. IVIII 
shows how to generate consistent approximations. The 
specialization to limiting cases, in particular to the dis- 
ordered but non-interacting electron system is given in 
Sec. IVIIII This conclude the general build-up of the the- 
ory. The case of disorder in the interaction part is briefly 
sketched in the Appendix [A] In the rest of the paper 
several concrete approximations are derived and classi- 
fied. This includes well-known but also new approxima- 
tion schemes. Mean-held approximations arc disussed in 
Sec. IIX[ cluster approximations in Sec. |Xj A summary 
and a discussion of general topics in Sec. IXII concludes 
the paper. 

II. HAMILTONIAN AND DYNAMIC 
QUANTITIES 

We consider a system of fcrmions in equilibrium at 
temperature T and chemical potential [i. In the grand- 
canonical ensemble the macrostatc of the system is given 
by the density operator 

cxp(-(g- M A)/T) 
9 tTexp(-(H-fxN)/T) ' 1 ' 

where N is the total particle-number operator, and H is 
the Hamiltonian. H is assumed to consist of a free (bilin- 
ear) part Hq which exhibits the (discrete) translational 
symmetries of an underlying D-dimensional lattice, a dis- 
order potential -ffdis; an d an interaction part H mt : 

H = H(t,T l ,U) = H (t)+H dis (v)+H- mt (U) . (2) 

The free part 

H a (t) = Y,t a pclc (3) 

a/3 

is characterized by a set of hopping parameters t a p where 
an index a labels the states of an orthonormal one- 
particle basis Typically, a refers to the sites x 
of the lattice as well as to some local degrees of freedom 
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(e.g. spin projection a =T,I), i-e. a = {x,a). The full 
hopping matrix with elements t a p is denoted by t. 
The interaction part 



H int (U) = - U a p Sl c^pOyCg 



(4) 



is a four-fermion point interaction and is specified by the 
(Coulomb) interaction parameters U a ps~ r The full set of 
interaction parameters is written as U for short. 
The disorder potential 



(5) 



is bilinear and given in terms of parameters rj a p which 
are random numbers with a joint probability distribution 
P(rj) with P(r}) > and / dr]P(ri) = 1. The configura- 
tional average for any quantity depending on rj is: 



(A) P = J d V P(v)A v 



(6) 



For the theoretical setup, Hdis(v) an< i H- m t(U) arc taken 
to be completely general. The construction of mean- 
field approximations will be most convenient for a local 
(Hubbard-type) interaction U and a diagonal (local) dis- 
order potential, r\ a ^ = 8 a f3rj a , with independent energies: 

Using the functional-integral formalism^ the grand 
potential, 

Qt,r),U = —T\nZ t ,r) : u , (7) 
and the partition function for a given configuration 77, 
Zt.r,,u = trcxp(-(iJ(t, rj,U) - fiN)/T) 

= J D^*exp(A tjT , iC r, ££ .) , (8) 

depend on the model parameters via the action 

n,<x{3 

-\ E U °f>*i f^drC^Qir^iT^siT) . (9) 

Here &,(iw n ) = T 1 / 2 J 1/T dr e^» T ^(r) (£(*w n ) = 

Tl/2 Io /T dT e ~ iu " T Ci T )) are Grassmann fields at the 
fcrmionic Matsubara frequencies iuj n = i(2n+l)7rT (with 
n = 0, ±1, ...). The configurational average of the grand 
potential is given by: 



Qt.p.u = J dr} P(ri)Q, t ,r,,u = (^t, V: u) 



(10) 



The subscript P indicates the dependence on the proba- 
bility distribution. 



For later purposes we need the free one-particle Green's 
function, 



G 



1 



t,o,o 



iuj n + fi — t 



(11) 



which is a matrix with the elements Gi,o,o,Q/3(* | ^'n)- The 
dependence of Gt.0,0 on t is indicated by the subscript. 
Dependencies on the chemical potential [i and the tem- 
perature T will not be indicated, fi and T are assumed 
to be fixed. Similarly, 



iu) n + fx — t — T) 



(12) 



denotes the free Green's function in the presence of 
the disorder potential. The action determines the full 
Green's function Gt^ v ,u the elements of which read: 

Gt,ri,U ,ap\iw n ) 

= ry [ D € D C€ct(iVn)£8(iv n )exp(A ttV! u !6e *) . 

^t,ri,U J 

(13) 

Finally, we introduce the (interaction) self-energy 

^t,n,v = G t Q — r) — G t v U = G t r} - G t n U . (14) 

^t,r),u depends on the configuration 77. The configura- 
tion independent (full) self-energy 

St,p,u = G^q - T~ P U (15) 

is defined with the help of the averaged Green's function 

T t ,p,u = (G t , v ,u)p ■ (16) 



III. CONFIGURATION-DEPENDENT 
SELF-ENERGY FUNCTIONAL 

The main idea of the self-energy-functional theory 
(SFT) is to express a thcrmodynamical potential as a 
functional of the (interaction) self-energy which is sta- 
tionary at the physical self-energy of the system. Vari- 
ation of the self-energy is achieved by taking trial self- 
energies from an (exactly solvable) reference system and 
varying its parameters. To be able to evaluate the self- 
energy functional (which in most cases is defined only 
formally), it is of crucial importance that the reference 
system shares with the original system the non-trivial 
part of the functional so that this can be eliminated. De- 
tails of the SFT are described in Refs. [HHSHi! In 
the following, we present a construction of the self-energy 
functional for a fixed configuration r). The construction 
is non-perturbative (i.e. does not refer to formal sums of 
skeleton diagrams) and allows for a generalization in the 
case of disorder (see Sec. IVII) . 
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To start with, we note that the action can be consid- 
ered as a functional of the (inverse) free Green's function 
(Eq. pD): 



n,a0 

~\ E U °P*i [^drCir^ir^ir^sir) . (17) 

Here G^ 1 is considered to be a free "variable" . The phys- 
ical action At . v .u,^* (Eq- ©) is obtained by evaluat- 
ing the functional Au.££* [■■■] at the physical inverse free 
Green's function G 1 = G^ 1 (Eq. (JI2J)), i.e.: 



(18) 



Note that a hat is used to distinguish functionals from 
physical quantities. Additional dependencies of a func- 
tional (parameters) are indicated by subscripts. 

In the same way, via Eq. (|7|) and Eq. |(HJ), the grand 
potential can be considered as a functional of G^" 1 , and 
one has: 



Qt,r),u = ^cr[G t)fJj0 ] , flt,o,u = &u[Gt,o,t 



(19) 



Again, one has to distinguish clearly e.g. between fit, 77,17, 
the exact grand potential of the model H (t, r], U) on the 
one hand, and 0[/[Gq ], a functional of the variable Gq 1 
on the other. The latter only acquires the value rit.77.f7 
if evaluated at Gq 1 = G^T 1 . 
The functional derivative 

1 Sh v [Gy 1 ] 1 SZu [Go 1 ] 

T SG^ ~ ZulG- 1 } 5G^ 

defines a functional Gu as 



= ~Gu [Go 1 ] 
(20) 



Gl7,ap [G 1 ] 
1 



Zu [Gq 1 ] 



[Go" 1 ] 



(21) 



which has the property 

Gu [Gt^fl] = Gt^.u , Gu [G^o, ] = G t ,o,u ■ (22) 

Namely, at the physical inverse free Green's function the 
functional integral in Eq. ([20)) defines the physical inter- 
acting Green's function, see Eq. (flU)) . 

Up to this point the derivations are standard. The 
decisive point in the construction of the self-energy func- 
tional is the following equation: 



Gu [G" 



£ = G 



(23) 



solved for G. This formal solution G = G[/[S] then de- 
fines a functional Gu [X] (which parametrically depends 
on U), i.e. we have 



Gu 



= Gum 



(24) 



for any S by construction. It is important to note that 
the functional Gf/[£] is universal, i.e. it does neither de- 
pend on t nor r/. If evaluated at the physical self-energy, 
the functional yields the exact Green's function 



Gu [S 



t,T),U 



= G 



t,T),U 



(25) 



since, by definition, Gu [St,rj,t/] solves Eq. (|24[) if 

Gu [£t,r,,u] 1 = G^ 1 - E t , 77,17- 

The final step is to use the (universal, i.e. t and 77 in- 
dependent) functionals Sljy [Gq -1 ] and Gt/[£] to express 
the grand potential as a functional of the self-energy. We 
define the functional as 



Fu\n 



n 



u 



Gum 



TrlnGf/p] (26) 



where TrA = TJ2 n E Q e l "» 0+ A aa {iuj n ). Eq. (U} and 
Eq. (2H) imply 



1 sF V m 

T 



= -Gt/[S 



(27) 



Hence, Gf/[S] is a "gradient" of the "scalar" self-energy 
functional Fu [S] . The physical meaning of Fu [£] is ob- 
vious when comparing with the original derivation of the 
SFT (cf. Ref. HI): Fu[S\ is the Legendre transform of 
the Luttinger-Ward functional^ 

Now, the grand potential can be considered as a func- 
tional of the self-energy: 



t.77,r7[ 



Trm- 



G 



Fu\n 



(28) 



t,77,0 



Two properties of this functional are very useful: 

First, at the exact self-energy X = St, 77,17 the self- 
energy functional yields the exact grand potential: 



^■t,r).u\^t,r),u] — ^t,77,L7 



since 



(29) 
(30) 
(31) 



from Eq. ([26j) and Eq. (J25J) . Hence, Eq. (J29J) follows from 
Eq. ((Til). 

Second, consider the derivative: 



^t,Tj,u[^t.ri.u] — Tr In Gt, v ,u + Fu[^t,tj,u] 
from Eq. ([25} and Eq. HID and 

Fu[-E w } = &U [Gt^o] - TrlnGt,77 lf 7 



The only purpose of this is to constitute a relation be- 
tween the variables G and S which may formally be 
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The equation 



G 



(33) 



is a (highly non-linear) conditional equation in the vari- 
able £ with parameters t, r/, U which is solved by the 
physical self-energy X = Ht.^.u- It is by no means 
straightforward to find a solution, however, since the 
functional Gt/[S] is not known explicitly but was con- 
structed in a formal way only. Obviously, this is equiv- 
alent to a search for the stationary point of the grand 
potential as a functional of the self-energy: 



ST, 



= 



(34) 



This establishes a very general variational principle with- 
out the need for an expansion in powers of the interaction 
strength, i.e. the construction is non-perturbative. 



IV. APPROXIMATIONS FOR A FIXED 
CONFIGURATION 

For the discussion of possible approximations, we first 
consider a fixed configuration 77. Then t + r\ is a fixed 
matrix of hopping parameters but without translational 
symmetries. The idea of the SFT is to construct ap- 
proximations by searching for the stationary point of the 
functional Eq. (|28[) on a restricted domain of trial self- 
energies. Trial self-energies are chosen from a reference 
system which shares with the original system the same 
interaction part Hi nt (U). In the Hamiltonian of the ref- 
erence system, 



H' = H (t') + H int (U) , 



(35) 



the bilinear part Ho(t') is varied arbitrarily. We have set 
rf = for the reference system. However, no transla- 
tional symmetry is assumed for the hopping t'. 

Since the interaction part is the same and due to the 
universality of the functional Fu[T], only the first term 
of the self-energy functional of the reference system, 



n t ',o,r/[E] =Trln 



1 



G t'fi,o ~ s 



(36) 



differs from the functional of the original system Eq. 
Combination of the functionals Eq. (|28[) and Eq. 
therefore gives 



t',0,U[ 



Tr In ■ 



1 



- s 



Trin 



r- 1 

"t',0,0 



(37) 



The not explicitly known functional i^jyp] cancels out. 



To search for the stationary point of the self-energy 
functional of the original system, we insert as trial self- 
energies the exact self-energies of the reference system: 
S = Sf.o.Lr- This yields a function of t', 



(38) 



Searching for the stationary point of flt^ v ,u(t') as a func- 
tion of t' means to search for the stationary point of the 
exact self-energy functional Eq. (f25)) on the restricted set 
of trial self-energies generated by the reference system 
with parameters t' . From Eq. (f29|) and Eq. (fT4)) for the 
reference system we have: 

ftt,r,,u(t') = £l#,o,u + Trln— I ^— 

G t . v .o - T>t>,o,u 

- TrlnGt',o )t /. (39) 

The important point is that the r.h.s. can be computed 
exactly if the reference system is an exactly solvable 
model. Specifying a certain reference system means to 
generate a particular approximation. Typically, a suit- 
able reference system can be found for lattice models by 
tiling the original lattice into clusters of finite size and 
by neglecting the inter-cluster hopping. 



V. STATISTICAL SFT 

To make contact with the statistical DMFT^ffl we 
consider a system with local interaction and local (and 
uncorrelated) disorder. The statistical DMFT treats the 
disorder part of the problem exactly while the (dynami- 
cal) mean-field approximation is used for the interaction 
part. Within the framework of the SFT, a mean-field ap- 
proximation is generated by a reference system in which 
all sites are decoupled. This implies that spatial correla- 
tions are neglected altogether in the computation of the 
self-energy The local (temporal) dynamics, however, can 
be optimized by introducing additional local degrees of 
freedom in the reference system. For the Hamiltonian of 
the reference system this means to introduce additional 
uncorrelated sites ("bath sites"), the on-site energies of 
which as well as their hybridizations with the original 
correlated sites are treated as variational parameters. 

To be explicit, the discussion is restricted to the 
Anderson-Hubbard model 



U^2n x1 n xL . (40) 



Here, x refers to the sites of a lattice, n xa = c xa c xai and 
r] x are independent random numbers distributed accord- 
ing to 



P{V) - \{p{ri x ) 



(41) 



with some density p(r/ x ). A reference system generating 
a mean-field approximation is H ' = Y^ x H x with (see 
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a) 



a) 



b) 



FIG. 1: (Color online) a) Representation of the Anderson- 
Hubbard model in one dimension. Squares/diamonds: fixed 
configuration of sites with local (Hubbard) interaction and 
local binary-alloy disorder. Lines: nearest-neighbor hopping, 
b) Reference system generating a mean-field approximation 
[n a = 4). Circles: bath sites (no interaction) with on-site en- 
ergies and hybridizations to original sites (square, diamonds) 
to be treated as independent variational parameters. 



Fig. mi 



(J X 

n s 

a t=2 



a i=2 



(42) 



It consists of effective impurity models with n s sites each: 
The correlated site x (with U ^ 0) and n s — 1 bath sites 
(with U — 0) labeled by i. The effective impurity models 
can be solved independently to get a trial self-energy. 

There is an indirect coupling, however, via the 
optimization of the variational parameters t' = 

(t' xx ,4 X) , V t X) y- The Euler equation, 



0_ 
di 



(43) 



simplifies due to the fact that the trial self-energy (and 
its derivative w.r.t t' , see Ref. 1571 ) is necessarily local. 
Using Eq. and carrying out the derivative, Eq. 
can be written as: 



1 



-Gv 



t'0,U 



n,x \ t,T),t 



dT, xx (iuj n ) 
dt' 



= 



(44) 



Now, the variation of the one-particle parameters of the 
impurity model at site x does not affect the self-energy 
of the impurity model at site x' (x' ^ x). Therefore, the 
Euler equation simplifies to: 




-Gv 



t'0,U 



dT, xx (iuj n ) 



= 0. 



(45) 



where t' xx denotes the variational parameters at the cor- 
related site x. This is a set of equations labeled by the 
site index x. Due to the matrix inversion in Eq. (|45[) . 



b) 



ooooooooo 



FIG. 2: (Color online) The same as Fig.[T]but for a continuum 
of bath sites (n s — > oo), represented by big circles. The ref- 
erence system generates the generalized or statistical DMFT 
(see text). 



however, the equations are coupled. This implies that 
for a generic configuration, the individual self-energies of 
the effective impurity models arc different at the station- 
ary point as the sites of the original lattice model are 
inequivalent. 

For n s — ¥ oo, i.e. for a continuous bath (see Fig. EJ, an 
additional formal simplification is possible: As it is obvi- 
ous from Eq. (|45|) . a solution t 1 of the coupled equations 



1 



Gt^fii^n) - St',o,t/(iw ra ) 



= (Gt',0,u(i^n))a 



_ (46) 

yields a stationary point of the self-energy functional. 
Note that for any finite n s < oo there is no solution: The 
r.h.s. of Eq. (|46|) is the Green's function of a finite system 
which exhibits a finite set of poles on the real frequency 
axis (after analytical continuation). Contrary, the l.h.s. 
represents an approximate lattice Green's function which 
in the thermodynamical limit L —> oo has branch cuts on 
the real axis induced by the branch cut of the free Green's 
function. 

In the case of systems with a few inequivalent sites, 
i.e. for inhomogeneous systems with a somewhat reduced 
translational symmetry, the equations (|46p exactly re- 
cover a generalization of the DMFT which has been put 
forward to describe correlation effects at surfaces and 
in thin films^ They are just the self-consistency equa- 
tions of this generalized DMFT. Typically, only a few 
effective impurity models have to be considered in this 
approachi 59 i 60 

For disordered systems without any translational sym- 
metry, the self-consistency equations f4"o]) constitute the 
statistical DMFT as introduced by Dobrosavljevic and 
Kotliar. 39 ! 40 In principle, the Eqs. (|46|) can be solved it- 
eratively. For any iteration in the self-consistency cycle, 
one then needs the local Green's function at each site 
of the system which in each case requires the solution 
of an interacting impurity problem. Eqs. (|46[) and the 
absence of translational symmetry also imply the need 
for an inversion of matrices with dimension given by the 
system size. For these reasons, the statistical DMFT is 
a numerically extremely expensive method. 

Choosing a reference system consisting of decoupled 
effective impurity models with finite (and actually small) 
n s could thus be an interesting alternative. Calcula- 
tions based on such a statistical dynamical impurity ap- 
proximation (stat-DIA), however, have not yet been per- 
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formed. Since all physical quantities derive from an ex- 
plicit though approximate expression for a thermody- 
namical potential, the stat-DIA is a thermodynamically 
consistent approximation. This can be seen as an ad- 
vantage compared to stat-DMFT approaches which em- 
ploy additional approximations to render practical calcu- 
lations possible. 

The self-consistency equations (|46p allow for a stochas- 
tic reinterpretation: For a given configuration of on-site 
energies {r/ x }, the local Green's function G xx will be site- 
dependent. The distribution of the local Green's function 
at a site x (generated by all configurations), however, will 
be the same as the distribution at a site x' ^ x since the 
on-site energies rj x have been assumed to be independent 
random numbers distributed according to the same den- 
sity p(r/ x ) for each x. Moreover, the distribution of the 
local Green's function at a site x, generated by all config- 
urations {rjx], is identical to the distribution of the local 
Green's function on all lattice sites for one fixed typical 
configuration of on-site energies. Hence, the index x in 
the self-consistency equations P5|) can be viewed as a 
label for a particular realization of the random variable 

( ' ;r. :i; • 

For a Bethe lattice, the equations (|46|) can be reinter- 
preted as stochastic recursion equations. Starting from 
an arbitrary initial sample for the local Green's function, 
{G ra }, the equations recursively generate a sequence of 
samples which converges to a sample which is represen- 
tative for the distribution of the local Green's function. 
The practical advantage of this reinterpretation consists 
in the fact that a Gibbs-sampling Monte-Carlo algorithm 
for the calculation of marginal distributions can be ap- 
plied (see Ref.|4l|, for an example). The iterative solution 
of the self-consistency equations for a given typi- 
cal configuration of on-site energies is then equivalent 
with the recursive update of a sample of local Green's 
functions. Furthermore, the matrix inversion required in 
Eq. pi>|) can be avoided in the case of a Bethe lattice. 
For a general (e.g. cubic) lattice, however, the stochastic 
scheme breaks down, and one has to return to the site 
interpretation again. 



VI. CONFIGUATION-INDEPENDENT 
SELF-ENERGY FUNCTIONAL 

The above discussion has shown the practical needs 
to construct more simple approximations. An intuitive 
strategy is to consider quantities involving configura- 
tional averages and to search for sensible approximations 
of the averaged quantities instead of considering full dis- 
tributions. The simplest and most natural average is the 
arithmetical average (•••)p which has been introduced 
in Sec. [Til It shall be understood that one has to be 
extremely careful when discussing transport properties 
in terms of the averaged one-particle Green's function 
^t,p,u : Close to Anderson localization the distribution 
of the local Green's function (at uj — 0) can exhibit an 



extreme asymmetry and a long tail such that the average 
is of no physical meaning and can by no means serve as 
an order parameter for a metal-insulator transition^ 

Here our goal is to construct non-perturbative and 
thermodynamically consistent approximations for aver- 
aged quantities which give information on thermody- 
namic properties and one-particle excitations. Within 
the framework of the self-energy-functional approach, 
this can be achieved by introducing functionals that in- 
volve quantities averaged according to the given probabil- 
ity distribution P. In particular, we consider functionals 
of the configuration-independent (full) self-energy S and 
the configuration-dependent (interaction) self-energy Yi^ 
as defined at the end of Sec. ITTI 

Analogous to Eq. (|23p . the equation 

( r- + s-„-s, ) p = r < 47 > 

constitutes a relation between the averaged Green's func- 
tion r on the one hand and S and H n on the other. The 

functional Qu [■ ■ ■] in Eq. ([2"3"]) is replaced by the func- 
tional (!/(■ ■ ■ — v))p> an d the probability distribution P, 
instead of the interaction parameters C7, plays the role 
of the external parameters. Contrary to Eq. ([2"3"j) . the 
above relation is diagonal in the frequency iu> n (which is 
suppressed in notations). 

Assume that S and (for any configuration rf) E,, are 
given. Then, the equation can formally be solved for V. 
This defines a functional Tp [■ ■ ■ , • • • ] which assigns the 
averaged Green's function T = Fp[S, to any S 

and S^. This functional plays a role analogous to the 
functional Gu [X] in Sec. IIII1 

Analogous to Eq. (|25p . we have: 

Fp{St,p,u, {^t,ri,u}] = ^t,p,u (48) 

since Eq. (jT7)) holds when evaluated for S = St.p,u, 
U, and r = Tt.p,u as it is obvious from Eqs. 
CE1]), (HI]) and HID. From Eq. §5$ we have 

GuP v ] = — 3-^ — (49) 

"■4,0,0 - "H ~ 

for any r\. This equation and 

f P [S,{S„}]= „_/_ g (50) 
"t.o.o 

form a (highly non-linear) system of conditional equa- 
tions for the variables S(iu) n ) and Ti^iiun)- The exter- 
nal parameters t, P, U specify the model under con- 
sideration. Eqs. (j49| and ([501 arc satisfied for the ex- 
act self-energies S(iu> n ) = St,p,u(i^n) an d S^(iw„) = 

St,r7,[/(*^n)- 

In the following we show that the conditional equations 
Eq. (J3U) and Eq. ([ST)]) can be considered as stationarity 
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conditions of the averaged grand potential as a functional 
of the self-energies. We define the self-energy functional 

1 



ftt,p,u[S, {^r,}] =Trln— T 



Tr In ■ 



r P [s,{n v }]-i + s-Ti--E Vi 



-Ttlnr P [S,{S ?7 }] + ( J F 1 t/ [S r ,] 



(51) 



The sum of the second and third term on the r.h.s. is 
a functional which is universal, i.e. it is independent of 
t (note that the terms do not cancel each other as the 
operations ln(- • • ) and (• ■ ■ )p do not commute). With 
Eqs. (0, |H]), Jill) and Eqs. (jUJ), (48]), the evaluation 
of the functional at the exact self-energies yields: 

Qt,p,u[St,p,u, {^t,ri,u}] = (TrlnGt, ?7 ^i7)p 

(52) 

i.e. the exact averaged grand potential. The functional 
derivatives are readily calculated: 



1 6nt,p,u[sdVr,}} 

T SS 
and 

1 Sht^ujS^-Er,}} 



1 



-r P [S,{V v }] (53) 



T 



5£ T 



Gu\E v ])p(ti). (54) 



Hence, setting the functional derivatives to zero, yields 
two equations equivalent with Eq. (|4T)|) and Eq. ([50")) . 
Therefore, the functional is stationary at the exact self- 
energies: 



SQt,p,u[St,p,u,{^t,r),u}] 

Is 

5Qt, P ,u[St. P ,u, {^t.j].u}] 



= 



. 



(55) 



The self-energy functional Eq. (|51[) represents a general- 
ization of the self-energy functional Eq. (|2"8"|) for interact- 
ing systems with disorder. It is completely general and 
provides an exact variational principle. 



VII. CONSISTENT APPROXIMATIONS 

In the spirit of the SFT for pure systems, approxima- 
tions shall be constructed by restricting the domain of 
self-energies in the functional Eq. (|51[) while retaining the 



exact functional dependence. We consider both, the full 
as well as the interaction self-energy. Trial self-energies 
are taken from a reference system which is a system in 
the same macroscopic state, i.e. with the same tempera- 
ture T and the same chemical potential /i as the original 
system, but has different one-particle parameters t' . The 
Hamiltonian of the reference system reads: 

H' = H(t',r),U) = JT (f ) + H dis (ri) + H lnt (U) . (56) 

H' has the same interaction part as compared to the 
original system. Likewise the disorder potential, i.e. the 
distribution P(r}), is assumed to be unchanged. Hence, 
the self-energy functional of the reference system is given 
by 



n tl , pu [S,{'E t1 }}=Trln 



Tr In ■ 



1 



Gt> ,o,o & 



-Trmr P [S,{S n }] + (j 1 [/ [S I7 ] 



(57) 



Only the first term on the r.h.s. is different as compared 
to the functional for the original system Eq. |5T|) . Com- 
bining Eq. ([ST]) and Eq. ([57]) . the last three terms on the 
respective r.h.s. cancel out, and one is left with 

1 



SltMlS^Vr,}} = JW^S,^}] + Trin — 

1 



Tr In ■ 



^t',o,o & 



(58) 



Note that the full and the interaction self-energies are 
considered as variables at this point, and that for the 
cancellation of the functionals it is of crucial importance 
to choose the reference system to have the same interac- 
tion and disorder. The self-energy functional Eq. ([53]) is 
still exact. 

As trial self-energies we insert the exact self-energies 
of the reference system: S = Sf , P .u an< i = ^t',rj,u- 
Searching for the stationary point of the exact self-energy 
functional Eq. ([5T|l on the subspace of trial self-energies 
taken from H' and parameterized by t' , means to search 
for the stationary point of a function of t': 

n t ,p,u(t') = &t,p,u[St>,p,u, {£tf,t|,u}] , (59) 

where t, P and U are fixed by the original system. From 
Eq. (J511) and Eq. (|15[) for the reference system we get the 
comparatively simple result: 



£lt,p,u(t') = Qf.p.u + Trln 1 



1 



G t,o,o _ Sf.p.u 
- Ti\sxT v ,p, v . (60) 

This result is formally very similar to Eq. (|39p for pure 
systems. Again, the important point is that the r.h.s. can 
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be computed exactly if the reference system is an exactly 
solvable model. The only difference consists in the fact 
that the grand potential, the Green's function and the 
self-energy of the reference system on the r.h.s. are re- 
placed by the corresponding averaged quantities and the 
configuration independent (full) self-energy. 

A certain approximation may be constructed along 
the following steps: (i) A reference system is specified 
with P and U fixed as in the original system. The 
hopping part, i.e. however, is fully at one's disposal 
and should be used to simplify the problem posed by 
the reference system. (ii) For a given set of varia- 
tional parameters t' , the reference system's Hamiltonian 
H' = Hit 1 ,T],U) is diagonalized for any configuration 
ry to get the (many-body) eigenenergies and eigenstates. 
(iii) The grand potential fit'. 77,17 and, from the Lehmann 
representation, the Green's function Gt\ v ,u are obtained 
for any rj. (iv) Averaging yields Qt',p,u = (&t',r},u)p 
and Tt'.p.u = (Gf ,v,u) p- The self-energy is computed 
via Sf.p.u = G^i q — p jj . (v) Inserting these results 
as well as the free Green's function of the original model 
into Eq. ([60]l yields Q t ,p,u{t')- (vi) Steps (ii) - (v) are 
repeated for different i! to find the stationary point t' s 
given by 



dn t ,p,u(iQ 

dt' 



= 



(61) 



This approximation strategy shares a number of advan- 
tageous features with the corresponding strategy for pure 
systems: Any approximation constructed in this way is 
a thermodynamically consistent one since the thcromdy- 
namics as well as the (averaged) one-particle excitation 
properties both derive from an explicit expression for the 
approximate averaged grand potential £lt,pu(t' s ) at the 
stationary point (see the discussion in Ref . \6u) . The only 
approximation consists in the restriction of the domain of 
the self-energy functional. The approach is systematic as 
an enlarged domain leads to an improved approximation 
(see the discussion in Ref. [57h . As the exact functional 
form is retained, approximations are non-perturbative by 
construction. 



VIII. LIMITING CASES 

Eq. (|5Tj) gives the self-energy functional 
Qt,p,u[Si {^r?}] for a disordered and interacting 
system. To discuss the limiting cases of the pure, of 
the non-interacting and of the pure and non-interacting 
system, we first have to specify the domain of this func- 
tional. This also applies to the functionals TpfS 1 , {Sj,}], 
G v m, etc 

For a given set of interaction parameters U and for 
a given probability distribution P, the domain T> = 
T>s x 2?£ of the self-energy functional Qt,p,u[S, {S,,}] 
shall consist of (full) self-energies S € T>s and (for any 
rf) interaction self-energies H v € T>s taken from the ref- 
erence system. Namely, a (full) self-energy S belongs to 



Vg, if there is some t' such that S = St>,p,u-, i-e. such 
that S is the exact self-energy of the problem given by 
H' = H(t' ,r),U) and P for some t' . Likewise, a set of 
interaction self-energies {£■,?} (for all possible 77) belongs 
to 2?s, if there is some i', such that X,, = Xf^u, i- c - 
such that Xjj is the exact (interaction) self-energy of the 
problem given by H' = H(t',r],U) for some t' . Hence, 
the hopping parameters t' span the domain of the self- 
energy functional. 

This definition is very convenient as it automatically 
ensures the correct analytical and causal properties for 
any self-energy in the domain. It also avoids formal dif- 
ficulties for pure or non-interacting systems which arise 
from the fact that conditional equations such as Eq. (|4"7|) 
become tautological in these limits and cannot serve to 
define a self-energy functional. With the above definition 
of the domain, however, this becomes irrelevant as for the 
cases of pure or non-interacting systems the domain con- 
sists of a single element or a null set only. 

With these preparations, let us discuss the limits in 
detail: 

(i) The pure and non-interacting case is given by 
P(f]) = 5(t]—t]o) and U = 0. Note that we can set rj = 
for simplicity (in the absence of disorder, a non-zero r/o 
merely implies a redefinition of the hopping: t + rjo). The 
domain of the self-energy functional Eq. (f5Tj) shrinks to 



the point S = S, 
(f25| . this implies 



0. According to Eq. and Eq. 



F u=0 [0] = h u=0 [G^lo] - Tr In G t . ,o = (62) 

as TrlnGt.0,0 is the grand potential of a system of non- 
interacting electrons with hopping t. For r\ = T70 = 
and S = H v = the second and the third term on the 



r.h.s. of Eq. (fSTj) cancel, and thus one is left with 
fit,fb,o[0,0] =TrlnG t , 0) o 



(63) 



which is the correct result. 

(ii) In the case of P(r]) = S(rj) but finite interaction 
U 7^ 0, one still has Y, v = S = S on the domain, and 
due to the cancellation of the second and the third term 
on the r.h.s. of Eq. (|5Tj) the self-energy functional reduces 
to Eq. (f2"5|) , as expected. 

(iii) For a system of non-interacting electrons (U = 
0) moving in a disorder potential with P{rj) 7^ S(r]), 
one has X,, = on the domain of the self-energy 
functional and thus F u=0 [0] = 0. With T P [S] = 
Tp[S, S,, = 0], this yields the self-energy functional 
fi^HS 1 ] = tot,p,u=o[S, ^ v = 0] with 

O t , P [S] = Trlii 1 - - -Trlnf P [S] 



Trln- 



T P [sy 



s-v, 



(64) 



for the problem with disorder only. The last two terms 
play the same role as the functional Fj/^P] for the prob- 
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lem with interaction only. This functional is discussed in 
the next section. 



IX. MEAN-FIELD APPROXIMATIONS 

Mean-field approximations for systems with local dis- 
order and local interactions represent a simple but in- 
structive class of approximations within the framework of 
the self-energy-functional approach. It is well known that 
any mean-field theory of disorder will be deficient in vari- 
ous ways. Issues such as localization cannot be addressed, 
for example, by means of the famous coherent-potential 
approximation (CPA).— Nevertheless, the mean-field 
concept represents an important benchmark and start- 
ing point for improvements and is in many cases the best 
we have at hand for practical calculations. 

We start by considering the functional (|64[) for the non- 
interacting, disorder-only limit of the model (|40|) . This 
is the Anderson model 

H = y^aa'C^Co:' + r] x cl s c a . (65) 

xx' x 

with local disorder given by Eq. (|4Tj) and some density 
p{r] x ) characterizing, for example, an alloy with R com- 
ponents (X) r Pr = !) : 

R 

Piv) = ^2Pr8{V - Vr) • (66) 
r=l 

For simplicity, the spin index is suppressed. 

A mean-field or single-site approximation is generated 
by a reference system consisting of decoupled sites (see 
Fig. [3b) , i.e. by switching off the hopping term. For an 
alloy-type disorder, this reference system is exactly solv- 
able, as one has to compute the Green's function for a 
finite number of R configurations. Models with a con- 
tinuous distribution p(r] x ) have to be simulated by a fi- 
nite but large R. Adding "bath" sites to the reference 
system, i.e. sites with fixed, configuration-independent 
on-site energies, enlarges the space of variational param- 
eters t! and trial self-energies S = St^p and implies an 
improved mean-field approximation (see Fig.[3b,d). Note 
that the disorder part is still the same as in the original 
model ([65]) . as it is required to justify Eq. (f60|) . 

The Hamiltonian of the reference system reads: 

H = ^ ] txx c x c x + T lx c x c x 

X X 

+ EE v t ] + h - c -) + E E ^"U- ■ 

x i x i 

(67) 

It consists of the local part of Eq. (|65[) and, for each site 
x, includes bath sites with (configuration- independent) 

(x) ... . . . . (x) 

energies e\ hybridizing with the original sites via 
where i = 2, ...,n s . As the Hamiltonian describes an 



b) 

OOOOOOOOO 

c) 5 5 5 5 5 5 J 5 5 
OOOOOOOOO 

FIG. 3: (Color online) a) Representation of the (disorder) 
Anderson model Eq. (|65[1 . b) Reference system generating 
the atomic approximation (see text), c) Two-site approxima- 
tion, d) Improved mean-field approximation due to atomic 
reference system with more bath sites (n s = 4). e) The opti- 
mum mean-field approximation generated by a continuum of 
bath sites is the coherent-potential approximation (CPA). If 
a) represents the (pure) Hubbard model, b) yields a Hubbard- 
I-type atomic approximation, d) is a typical mean-field ap- 
proximation, and the optimum mean-field approximation e) is 
given by the DMFT. In case that a) represents the Anderson- 
Hubbard model Eq. (|40p , the reference system d) leads to the 
DMFT+CPA approach. 



impurity model with identical and decoupled replicas at 
any site a;, it is in fact sufficient to focus on one impurity 
model only. The site index x can be suppressed in this 
case. This reflects the translational symmetry of aver- 
aged quantities in the original model (|65|) . 

Due to the decoupling of the original sites, the ref- 
erence system yields a trial self-energy which is local. 
Spatial correlations due to non-local contributions of the 
self-energy are neglected. Differences between different 
mean-field approximations are due to the temporal cor- 
relations, i.e. due to additional bath sites. Obviously, 
the optimum single-site approximation is obtained for a 
continuum of bath sites n s — > oo (Fig. [3k). 



A. Atomic approximation 

The simplest approximation is obtained for n a = 1, 
i.e. no additional bath sites (Fig. (3b) ■ This case is in- 
structive as it allows for a largly analytical treatment 
which elucidates some general features of the disorder 
SFT. The reference-system Hamiltonian consists of the 
first two terms on the r.h.s. of Eq. ([67]) . We consider 
the case of a binary alloy (R = 2) with on-site ener- 
gies rji and 772 and corresponding probabilities p\ and pi. 
The only variational parameter left is the (configuration- 
independent) on-site energy t' = t' xx . 

For T = and /i = the averaged grand potential of 
the reference system tt' = flf,p,o is easily calculated: 

si'/l = Pi{t' a + m)Q(-t'o - vi)+P2(t' + V2)Q(-t' - m) ■ 

(68) 
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This is the first term on the r.h.s. of Eq. ([60)) . 

The averaged Green's function of the reference sys- 
tem (Eq. (fT6|l ) is obtained immediately as T' mas ,(w) = 
8 xx 'V(lj) with 



r» = 



pi 



P2 



W — t' Q — Tjl LU — t' Q — 7/2 



(69) 



Using the considerations of Ref.[27| (see Eq. (20) therein), 
only the poles uj^ 2 = t' +r)i 2 enter the result for the third 
term on the r.h.s. of Eq. (J5D]): 



TrmIY,P,o/L= J2 v'M-w't)-Rs/L. 



(70) 



The term Rs (see Ref. |27|) cancels out later. 

With the free Green's function of the reference system, 
Gq(uj) = l/(ui — t' ) we get from Eq. (fT5|) the self-energy 
as S xx '(lo) = S xx 'S(uj) where 



{if) - {Tff 



t'a - m - m + (v) 



(71) 



Let e(fc) be the eigenvalues of the hopping matrix t and 
let T xx ,(uj) = L- 1 J2 k e tk{x - x,) T k (uj) with L the number 
of lattice sites and Tk(uj) = 1/{uj — e(k) — S(cj)) be the 
averaged Green's function of the original system as it 
appears in Eq. (|60[) . From Eq. (|7Tj) we find 



Oil 



a 2 



UJ2 



(72) 



with poles wi, 2 = (e(fe) +t' + 771 +r? 2 )/2± [(e(fe) + 2(77) - 
Vi <o) 2 / 4 + (V 2 ) - (v) 2 } 1/2 and weights a x = (w x - 

?7i — 772 + <?7> — *o)/ (^1 —^2) and a 2 = (uj 2 - Vi - m + (v) - 
t' )/(u>2 —uii). The k dependence of the poles is only due 
to s(k), i.e. uv = Lj r (k) = u> r (e(k)). Hence, using Eq. 
(21) of Ref. |23 and the definition of the free density of 
states po(z) = L^ 1 J2 k 8{z — e(fe)), we get 

/oo 
dz p (z)aj r (z)@(-uj r (z)) 
-OO 



r=l,2 



- Rs/L. 



(73) 



for the second term on the r.h.s. of Eq. (l60l) . 

Adding the three contributions, one can search for 
a stationary point numerically. Here, we restrict our- 
selves to the particle-hole symmetric case, i.e. we assume 
p (z) = Po(~z), m = m and p 1 = p 2 = 1/2. While 
the different contributions to the SFT grand potential 
Eq. (pO)) are asymmetric, it is straightforward to see that 
their sum is symmetric with respect to a sign change of 
the variational parameter: ^lt,p,o{t'o) = ^t,p,o{~t'o)- Fur- 
thermore, the dependence on t' is smooth for \t' \ < \rji\. 
Hence, the SFT grand potential is stationary at t' Q = 
which had to be expected by virtue of particle-hole sym- 
metry. 

The optimum self-energy is given by Eq. (|7ip with 



t'n 



0. It consists of the "virtual-crystal" potential 



(77) (note that S(u>) = (77) is the so-called virtual-crystal 
approximationii) and a frequency-dependent part with 
one simple pole at w = 0. Its weight is given by the dis- 
order strength, namely by the variance (77 2 ) — (?7) 2 = r/f. 
For any finite 771 this leads to a splitting of the aver- 
aged local density of states into a lower and an upper 
alloy band. In the case of strong disorder, this result is 
qualitatively correct as could be expected since the high- 
frequency behavior of S(u>) is correct up to the order 
1/ui 2 . However, the widths of the alloy bands turns out 
to be too small. 

This atomic-like approximation is very much remi- 
niscent of the Hubbard-I approximation^ for the pure 
but interacting system - although the Hubbard-I self- 
consistency condition is somewhat different from the SFT 
Euler equation and leads to different results away from 
the particle- hole symmetric point. The analogy between 
approximations for disordered but non-interacting and 
pure but interacting systems relies on the same struc- 
ture of the reference system. In this analogy, the static 
part of the disorder self-energy, the virtual-crystal poten- 
tial (77), corresponds to the static part of the Hubbard- 
I self-energy which, for the Hubbard model, is given 
by U(c xa .c xrJ ), i.e. to the Hartree(-Fock) approximation. 
A combination of both Hubbard-I-typc approximations 
for the interacting and disordered system is straightfor- 
wardly set up with an atomic-like reference system in- 
cluding (local) interaction and disorder. 

In view of the simplicity of the Hubbard-I-type ap- 
proach, it is remarkable that the variational optimization 
of the on-site hopping ensures thcrmodynamical consis- 
tency with respect to the particle number, i.e. the av- 
eraged particle number as obtained from the (approxi- 
mate) disorder-averaged grand potential as a p, deriva- 
tive is always the same as the averaged particle number 
calculated by integration of the (approximate) disorder- 
averaged single-particle Green's function. The proof for 
this consistency is analogous to that given in Ref. [6l|, see 
Appendix [B] 

It is also interesting that the stationary point t' = 
actually (locally) maximizes the SFT grand potential as 
has been verified by a simple numerical evaluation of Eq. 
([6H)l for the particle-hole symmetric case. In general, and 
for high-dimensional parameter spaces in particular, we 
expect that stationary points are saddles. 



B. Two-site approximation 

The two-site approximation (see Fig. [5J:) is the sim- 
plest mean-field approach beyond the atomic approxi- 
mation. For the case of the pure Hubbard model with- 
out disorder is has proven to be very instructive and 
successfuh 27 i 63 ' 64 i 65 ' 66 It provides a handy mean-field 
approach which is able to reproduce qualitatively the 
DMFT phase diagram for the Mott transition in the 
single-band model and which has been employed to study 
more complex two- and multi-orbital systems. A recent 
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application^ to the Anderson-Hubbard model has shown 
that the approach can straightforwardly be extended to 
disordered (and interacting) systems along the lines de- 
scribed here. The two-site approximation qualitatively 
reproduces the results of the DMFT+CPA method but 
with a minimum computational effort. 



C. Coherent-potential approximation 

The two-site approach can of course be improved by 
adding more bath sites (Fig. [3ji). The best mean-field 
approach is then obtained for n s — > oo, i.e. for a contin- 
uum of bath sites (Fig. (3^)- Varying all parameters, i.e. 



t' = t' 



ei x) and Vi 



V: 



(as) 



yields the optimum 



local self-energy as the stationary point 



n d ri rc i V^V^ 3^t,p[St',p] dS xx (iuj r , 

= = 



8S xx {iui n ) 



at' 



(74) 



of the function Q, t p(t') = f2t p[St' p]- Here it has been 
used that the self-energy (and also its i'-derivative) is 
local and non-zero at the impurity site only. From Eq. 
([6"4"| we have: 



n t ,p(t') 



Tr In 



1 



G t,o,o ~ s *',p 



Tr \nT P [St>,p] 



Tr In ■ 



1 



T P [S t ,, P }-^ + S t ^ P -r 1i 



.(75) 



Note that the last term on the r.h.s. is just the averaged 
grand potential of the reference system flt^p. 
The Euler equation (|74j) is satisfied, if 







8Q t ,p{St> ,p] 
SS xx (iu n ) 



(76) 



for each site x and for each Matsubara frequency u> n . 
Calculating the derivative of the functional Eq. (|6"4")h we 
obtain: 



1 



1 Sflt,p,o[St' ,p] 
T SS xx (iuj n ) \iw n + fx — t — S(iu) n ) 

(I 1 



T P [St'.p}- 1 + S(iuj n )-r li 



(77) 



The first term on the r.h.s. is the local element of the 
averaged Green's function of the lattice model, Ti oc (iu> n ) 
which is calculated with the approximate local self-energy 
S{iw n ) (exploiting translational symmetry, the site index 
can be supressed). The second term on the r.h.s. is the 
averaged Green's function of the reference system at the 
impurity site T[ oc (iui n ). The optimum local disorder self- 
energy is thus determined by the condition that the local 



averaged Green's function equals the averaged impurity 
Green's function of the reference system: 



T\ oc {iu n ) = r( oc (iu;„) 



(78) 



This is exactly the self-consistency condition of the 
coherent-potential approximation. 

To elucidate this point, we note that, for the reference 
system Eq. ([57)1 . T[ oc (iu) n ) is the average of the impurity 
Green's function for the different on-site energies, i.e.: 



r' 



loci 



= &r 



1 



■ fl — Tj r — A(iLU r] 



(79) 



where A(iuj n ) = ^ V? / (iiv n + fi— Ei) is the hybridization 
function. On the other hand, from the definition of the 
disorder self-energy for the reference system, we have 



r ioc(* w ») = 



1 



iuj n + n — A(iuj n ) — S(iu> n ) 
Eliminating the hybridization function, we get: 

1 



TlocOn) = 



Pr 



(80) 



(81) 



After a few manipulations, this equation can be cast into 
the form 



r) r - S(iui n ) 



T[ oc (iu) n )(r) r - S(iu) n )) 



(82) 



which makes contact with the original derivation of the 
CPA where the averaged atomic scattering matrix is set 
to zero^ Introducing the free Bloch-band dispersion 
e(fe) as the Fourier transform of the hopping t and the 
free Bloch-density of states po(z) = {^1 L)^ k 5{z — e{k)), 
we have 



rioc(«w n ) 



po{z)dz 



+ fi — z — S(iuj n ) 



(83) 



In combination with Eq. (J78J), this equation can be used 
to eliminate T[ oc (iu! n ) from Eq. (|8"Tj) to obtain a single 
conditional equation for S(iu! n ). 

For a practical calculation, one may set up the fol- 
lowing self-consistency cycle: Starting from a guess for 
the hybridization function A(iu) n ), the averaged impu- 
rity Green's function T[ oc (iu> n ) can be computed from Eq. 
([75]). With the help of Eq. flSC]) this determines S(iw n ) 
which is then used in Eq. (|83|) to get the CPA Green's 
function T\ oc (iu) n ) = T[ oc (iu> n ). Using Eq. (|5D|) again, a 
new hybridization function can be found. 

Obviously, the bath sites of the reference system play 
the role of an "effective medium". The bath parame- 
ters or, equivalently, the hybridization function A(iui n ) 
parameterize the local disorder self-energy in the most 
general way consistent with causality requirements. The 
present rederivation of the CPA therefore very clearly 
shows the CPA to be the best local approximation. 
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D. Dynamical mean-field theory 

The reference systems shown in Fig. [3] can also be used 
in the context of a pure system without disorder but with 
local interaction such as the Hubbard model. In this case 
the theory reduces to the conventional SFT. The refer- 
ence system b) with the hopping between the correlated 
sites switched off generates an atomic approximation very 
much the same as the Hubbard-I approximation^ but 
with the Hubbard-I self-consistency replaced by the SFT 
Euler equation for that reference system which is differ- 
ent. The rather crude atomic- like approximation can be 
improved by adding uncorrelatcd bath sites. This yields 
the reference system d) which is a set of disconnected 
single-impurity Anderson models with n s < oo sites 
each. Qualitatively, the results of the two-site dynamical- 
impurity approximatio n 27 ' 63 i 64 ' 65 i 66 (Fig. [3}:) are already 
close to the n s = oo limit. The convergence with increas- 
ing n s is rapid (see Ref. [(38|, for example). 

The most general local trial self-energy compatible 
with the requirements of causality is generated in the 
limit n s = oo. Here, the calculation proceeds in a way 
analogous to the previous section and yields the Euler 
(or self-consistency) equation (see Ref. \2w : 



Gl oc {iLU n ) = G[ oc {iuin) 



(84) 



Here Gi oc (iu n ) is the on-site element of the lattice 
Green's function (Fig. [3k,) calculated approximately from 
the approximate local self-energy via the Dyson equation 
for the lattice model, and G[ oc (iui n ) is the exact Green's 
function of the reference system (Fig.[3k) at the impurity 
site. Eq. (|84|) is just the self-consistency equation of the 
dynamical mean-field theor y 12 i 14 i 15 i 16 

The CPA for a non-interacting system with local dis- 
order has a formal structure which is very similar to the 
DMFT for a pure system with local interaction. This is 
apparent when comparing the respective self-consistency 
conditions (|78[) and (j84"|) and also the respective self- 
consistency cycles which serve to iteratively solve the 
mean-field equations. As the DMFT, the CPA becomes 
exact in the limit of high spatial dimensions D — > oo as 
has been shown by Vlaming and Vollhardti 13 ' 69 Both ap- 
proaches are characterized as approximations that yield 
the optimum local self-energy. All this is made very ob- 
vious within the framework of the SFT which discloses 
the formal analogies between Green's-function-based ap- 
proaches to disordered or interacting systems. 

A self-evident idea that has been persued in the past 
is to derive the DMFT (or a different many-body ap- 
proach to pure interacting systems) by using the formal 
structure of the CPA. This requires, however, a trans- 
formation of a pure system with local interaction to a 
non-interacting system with local disorder to which the 
CPA can be applied. Hubbard's alloy analogy 7 -^ and also 
a refined version 7 -^ represent such transformations. The 
subsequent application of the CPA to the Hubbard's fic- 
ticious alloy yields the so-called Hubbard-III approxima- 
tion. The alloy analogy (the transformation) itself, how- 



ever, must be seen as a rough approximation which fails 
to recover Fermi-liquid properties even for weak interac- 
tions. The "many-body CPA" by Hirooka and Shimizu, 72 
the "generalized CPA" by Janis 7 - 3 - as well as the "dy- 
namical CPA" of Kakehashi 7 - 4 - go beyond a simple anal- 
ogy. The idea of the dynamical CPA is to perform the 
transformation to an effective one-particle Hamiltonian 
within the functional-integral formalism by means of a 
Hubbard-Stratonovich transformation and to recover the 
DMFT by the subsequent use of the CPA (see Ref. [ll for 
a discussion). 

In this context it is also worth mentioning the Falicov- 
Kimball model which can be considered as a variant 
of the Hubbard model with the hopping of one of the 
two spin species switched off. The exact solution of the 
Falicov-Kimball model in the limit D ■— > oo has been 
worked out by Brandt and Mielschi 76 ' 77 i 78 They could 
show that the local (interaction) self-energy of the mo- 
bile carriers is given by the CPA (disorder) self-energy 
- no alloy analogy is necessary for this simplified model. 
It is tempting to understand the dynamics of the mo- 
bile carriers as the scattering of non-interacting particles 
from the (local and uncorrelated binary) disorder poten- 
tial generated by the immobile ones. 



E. DMFT+CPA 

For a system with Hubbard-type interactions and local 
disorder, e.g. the prototypical Anderson-Hubbard model 
Eq. (|40p , the optimum mean-field theory is generated by 
the reference system shown in Fig. [3^. Note that the 
reference system shares with the original system the same 
interaction part and the same disorder potential. The 
continuum of bath sites is uncorrelatcd and configuration 
independent. 

We start from Eq. ([60)) . For the present case the Euler 
equation dilt.p.u(t')/dt' = is satisfied if 



1 S£lt,p,u[St',p,u, {^t',Tj,u}] 



T 



5S xx (iu) n ) 







(85) 



for the local elements of the self-energy. This is the Euler 
equation which fixes the variational paramters t' . Anal- 
ogous to Eq. (|77|) we then obtain 



Tp[St',P,U,{^t',r),u}] 3 



fJt — t—S; 



t'.P.U 



(86) 

i.e. the local averaged Green's function equals the aver- 
aged impurity Green's function of the reference system, 
ri oc (iu;„) = r[ oc (ia> n ), as in the non-interacting case with 
disorder. 

Analogous to the procedure described in Sec. lIXCl this 
self-consistency equation can be solved in an iterative 
manner: Starting with a guess for the variational param- 
eters t' or, equivalently, for the hybridization function 
A(iw n ) = V?/(iui n + /j, — Si), the interacting impurity 
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Green's function is calculated for any (local) configura- 
tion of the reference single-impurity Anderson model and 
averaged over the configurations to get T[ oc (iuj n ). The 
self-energy is obtained from Eq. (|5D|) which is the defin- 
ing equation Eq. (| 1 5[) for the self-energy also in the case 
of an interacting impurity model. S(icu n ) is then used 
in Eq. to get the averaged lattice Green's function 
T\ oc (iuj n ). Via the self-consistency equation, this gives us 
r[ oc (iu) n ) and, using Eq. ([80|) again, a new hybridization 
function. Note that for a continuum of bath sites the op- 
timum on-site element of the hopping t' is always given 
by t'o = txx- This can be shown in essentially the same 
way as has been done in Ref. [6l| for the pure case. 

The self-consistency condition Eq. (f5o| and the cy- 
cle desribed above constitutes what is known as the 
DMFT+CPA approach which has been put forward 
by Janis and Vollhardt 3 ^ and by Dobrosavljevic and 
Kotlia r 31 i 32 in a different context: Originally, the 
DMFT+CPA has been introduced and charactized as 
the exact theory in the limit of infinite spatial dimensions 
which remains non-trivial for a proper scaling of the hop- 
ping parameters^ The presented rederivation places the 
DMFT+CPA into the broader framework of the SFT. 

The typical medium theory (TMT), as suggested re- 
cently by Dobrosavljevic, Pastor and Nikolic,— can be 
seen as a variant of the DMFT+CPA (see also Refs. 
I48ll49h . Here a modification of the above-described self- 
consistency cycle is considered by replacing {ad hoc) the 
usual arithmetical average of the interacting impurity 
Green's function over the (local) configurations by a geo- 
metrical one. Actually, this geometrical average is ap- 
plied to the corresponding (positive definite) spectral 
density, and the Green's function T[ oc (ioj n ) is obtained 
afterwards from this average as the usual Hilbert trans- 
form. A rederivation of this variant within the SFT 
does not seem to be possible: Any approximation gener- 
ated within the SFT preserves the elementary sum rule 
/ duo A\ oc {uj) = 1 for the local spectral density while this 
must be violated when averaging geometrically. Despite 
this conceptual shortco ming , the TMT can be motivated 
physically (see Refs. I47ll48[ ) and clearly improves upon 
the mean-field concept as it is able to describe aspects 
of Anderson localization. It should be seen as a prag- 
matic simplification of the statistical DMFT that allows 
for practical computations in many cases which are not 
accessible to the statistical DMFT. As can be shown for 
a simplified model^ the results of TMT qualitatively 
agree very well with those obtained from a full evaluation 
of the statistical DMFT for interaction- and disorder- 
driven metal- insulator transitions. 



X. CLUSTER EXTENSIONS 

A mean-field approximation neglects spatial correla- 
tions by considering an effective single-site problem to 
generate the self-energy. A straightforward idea to im- 
prove upon the mean-field concept is therefore to re- 
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FIG. 4: (Color online) Decomposition of real-space lattice 
vectors, x = x + X , and reciprocal-space wave vectors, k = 
k + K , for a D = 1 dimensional lattice (lattice constant a) 
with L = 12 sites tiled with L/L c = 3 clusters consisting of 
L c — 4 sites each, x: original lattice, x: superlattice. X: 
sites in a cluster. Reciprocal space: There are L allowed wave 
vectors k in the unit cell of the lattice reciprocal to x, and 
there are L/L c allowed wave vectors k in the unit cell of the 
lattice reciprocal to the superlattice x. K are the reciprocal 
superlattice vectors, exp(iK x) = 1. 



place the effective single-site model by an effective cluster 
model consisting of a finite (small) number of sites. This 
should restore the important short-range correlations and 
provide a systematic approach with the inverse cluster 
size as a small parameter. The goal is to construct a 
cluster approximation which on the one hand reduces to 
the CPA or to a simpler mean-field approach, which on 
the other hand becomes exact in the infinite-cluster limit 
and which preserves the translational and point-group 
symmetries of the lattice in addition. This, however, is 
by no means trivial, and naive theories often suffer from 
causality violations*^ In contrast, it is easy to see that 
all approximations that are constructed within the SFT, 
including the different cluster approximations discussed 
in the following, are manifestly causal. This is briefly 
discussed in Appendix [Cl 

Before discussing the different cluster approaches in 
detail, let us introduce some notations (see Fig. [4]): We 
consider a system on a D-dimcnsional lattice of L sites 
with periodic boundary conditions and L — > oo in the 
end. The position vector to a site in the lattice is denoted 
by x. There are L allowed wave vectors in a unit cell of 
the reciprocal lattice which are denoted by k. The lattice 
is tiled with L/L c clusters consisting of L c sites each. 
Let x be the position vector of the cluster origin, and X 
the position vector of a site in a cluster, referring to the 
cluster origin. We then have the unique decomposition 
x = x + X. The vectors x form a superlattice with 
a unit-cell volume enlarged by the factor L c . In a unit 
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cell of the reciprocal superlattice there are L/L c allowed 
wave vectors k. Its volume is reduced by the factor L c 
as compared to the volume of the reciprocal unit cell of 
the original lattice. For a given k we have the unique 
decomposition k = k + K where K are the vectors of 
the reciprocal superlattice, i.e. exp(iKx) = 1. In the 
reciprocal unit cell of the original lattice, there are L c 
vectors K . These can also be interpreted as the allowed 
cluster wave vectors when imposing periodic boundary 
conditions on the individual cluster. 

Consider the L x L matrix U with elements 



1 



(87) 



and the L/L c x L/L c matrix V with elements 



and the L c x L c matrix W with elements 

1 



W 



X.K 



iKX 



(89) 



U, V and W arc unitary and define Fourier trans- 
formations between the respective real and reciprocal 
spaces. Note that U ^ VW = WV. A quan- 
tity A XtX i which is invariant under lattice translations 
x , i.e. A x+X(hX , +Xa = Aa,^', is diagonalized by U: 
(W AU)kk' = A(k)5k,k'- A quantity A X:X > which is 
invariant under superlattice translations Xq as well as 
under cluster translations Xq (i.e. which is cyclic on 
the cluster), A x+Xo . x , +X() = A x+ x Q , x '+x a = A x ^ x >, 
is diagonalized by VW: (WW* AVW) %K £, K , = 

A{k,K)5 l% ,5 K , K <- 



A. Variational cluster approach 

A straightforward extension of the single-site atomic 
approximation (see Sec. IIX A[) is the (disorder) varia- 
tional cluster approach (VCA) which is the analog of 
the VCA know n 28 ' 53 ' 54 for the interacting but pure sys- 
tem. To be definite, consider the Anderson model Eq. 
([55]) again. The model is represented in Fig. [5h.. A suit- 
able reference system to include short-range correlations 
in the disorder self-energy consists of a set of isolated 
clusters of L c sites tiling the original lattice as shown by 
Fig. 03. Since within the SFT, the (only) approximation 
is to replace the exact self-energy by the self-energy of 
the reference system, the VCA must become exact in the 
limit L c — > oo, i.e. the VCA is a systematic approach 
which is controlled by the inverse cluster size as a small 
parameter. 

To calculate the cluster self-energy St\p, the Ander- 
son model must be solved on an individual cluster (by 
switching off the inter-cluster hopping) for all configura- 
tions. For large cluster size L c , the computational cost 
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FIG. 5: (Color online) a) Representation of the Anderson 
model Eq. ((65} . b) Reference system generating the (disorder) 
variational-cluster approximation (VCA) (without variational 
optimization: cluster-perturbation theory), c) Reference sys- 
tem generating the molecular CPA (M-CPA). In the case of 
an additional local (Hubbard) interaction, a) represents the 
Anderson-Hubbard model, b) yields a generalized VCA, and 
c) generates the C-DMFT+M-CPA approximation. 



therefore scales exponentially with L c . This is, of course, 
characteristic for all cluster approximations. 

The intra-cluster hopping parameters t' including the 
on-site energies are considered as variational parame- 
ters. A non-variational variant of the VCA is obtained if 
one sets these parameters to the original (physical) hop- 
ping parameters within the cluster: t' = tmtra- This 
constitues the analog of the cluster-perturbation theory 
(CPT) 51 i 52 which is known for the interacting, pure sys- 
tem. 

Obviously, the CPT is also systematic and controlled 
by l/L c . The additional optimization of the variational 
parameters t' within the VCA is expected to speed up 
the convergence with increasing L c as it does for inter- 
acting but pure systems^ One should note, however, 
that the variational optimization should become less and 
less important with increasing L c as the exact solution 
is approached for L c — > oo with t' = ti n t ra anyway. It 
is therefore advisable to restrict the space of variational 
parameters and to optimize a few parameters only as, for 
example, those parameters that are close to the cluster 
boundary. 

Parameter optimization is beneficial if (local) interac- 
tions are considered in addition (as e.g. in the Anderson- 
Hubbard model). Fig. [SJd defines a generalized VCA for 
this case. The additional interactions can drive a spon- 
taneous breaking of a continuous symmetry. This is sig- 
naled within the VCA by a stationary point of the SFT 
grand potential at a non- vanishing value of a symmetry- 
breaking field. Note that the respective field term can 
be added to the reference-system Hamiltonian if this is 
given by a one-particle term. This "Weiss" field is a fic- 
ticious one which clearly has to be distinguished from a 
physical field and which describes spontaneous opposed 
to induced symmetry breaking. For pure systems and 
spontaneous SU(2) and U(l) symmetry breaking, this 
concept has already been applied successfully. 53 ' 54 ' 61 ' 80 
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B. Molecular CPA 

The (disorder) VCA is obtained as the cluster general- 
ization of the atomic approximation. Likewise, within the 
framework of the SFT, the cluster generalization of the 
CPA leads to the so-called molecular CPA (M-CPA) 41 
To be definite, we again consider the Anderson model 
Eq. (|65|) . The reference system generating the M-CPA is 
shown in Fig. [SJ:. It consists of a set of isolated clusters 
of L c sites each and a continuum of bath sites attached 
to each of the original sites. 

As in the case of the CPA, this allows to derive a sim- 
plified Euler equation. Analogous to Eq. (|7^|) and Eq. 
(|77[) we can derive the following stationarity condition: 



= r EE 

n x,x' 

- T t ',p{iu> n ) 
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' \iuj n + /i — t — S(iuj n ) 
dS xx >(iu n ) 



x' X.ll 
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(90) 



x and x' must belong to the same cluster since 
S xx '(iuJ n ) = and also the "projector" dS xx '(iuj n )/dt' = 

if a;, a;' belong to different clusters as these are decou- 
pled in the reference system. Note that t' is a matrix 
labeled as t' xx , u , where i = 1, ...,n s and i ^ 1 refers to 
the additional bath sites attached to each original site 
x. The same holds for IV. p and S = St'.p- As the 
reference system exhibits local disorder on the original 
sites x, there are non-zero elements of the self-energy for 

1 = i' = 1 only. We write S xx '_\i = S xx i for short. 

Obviously, S xx > = = S X x>5 xx >. 

Both, the original and the reference system, are invari- 
ant under translations of the superlatticc. Hence, Fourier 
transformation given by V~ r (Eq. (|88[)) is appropriate. 
This yields: 
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(iLU n ) 



X'X.ll 



ioj n + fx — t(k) - S(iuj n ) 

dSxx'{iu n ) 
dt 1 



Here, k runs over the L/L c wave vectors in the reduced 
Brillouin zone. The trial self-energy S as well as IV, p 
are k independent matrices in the intra-cluster position 
vectors X,X' while t(k) is k dependent. 

The M-CPA self-consistency equation reads^i 
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where the frequency dependence has been suppressed for 
convenience. This generalizes the CPA self-consistency 
equation ([75]) which is obtained from Eq. (f9"2"|) for X = 
X' and cluster size L c = 1. 

Comparing with the stationarity condition Eq. (|9ip , we 
note that the SFT grand potential is stationary at the M- 
CPA sclf-cncrgy, i.e. at the self-energy S(iuj n ) = Sf,p if 
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FIG. 6: (Color online) Systematics of different approxima- 
tions that can be constructed within the SFT framework for 
pure interacting systems, disordered free systems or interact- 
ing systems with disorder. The interaction and/or disorder is 
assumed to be local. The axis intercepts L c and n s character- 
ize the reference system and thereby the type of approxima- 
tion. A local approximation ("dynamical impurity approx- 
imation", DIA) is obtained for L c = 1. The optimum local 
approximation with n 3 — > oo bath sites is the DMFT, CPA, or 
DMFT+CPA, respectively. A cluster ("cellular" or "molecu- 
lar") extension of DMFT/CPA is obtained for L c > 1. The 
variational cluster approximation (VCA) is a cluster approx- 
imation without the use of bath sites (L c > 1, n 3 — 1). The 
exact solution would be obtained in the limit L c — » <x irre- 
spective of n a . 



variational parameters t' can be found such that Eq. (|92|) 
is satisfied for any iuj n , X and X'. Hence, the reference 
system Fig. [SJ; generates the M-CPA. 

The following self-consistency cycle can be set up: We 
start with a guess for the one-particle parameters of the 
reference system, i.e. the parameters for a single cluster 
t' 



0,X,X' 



t' S+ x, S+X ',xi, 4 X) and F/ X) (i = 2,...,n s ). 
For a fixed (intra-cluster) configuration of the disor- 
der potential r]x, the intra-cluster part of the Green's 
(91) function Gt' , v is an L c x L c matrix G with elements 
By solving its equation of motion, one easily 



G 



x+X,x+X' 

verifies G'(iLu n ) = 1 / (icjn + /i — 1' — rj A.(iuj n )) where the 
matrix hybridization function A(icj„) is diagonal with 

elements A x {iu n ) = X^2 n I ( iuj n ^- X) )- Av- 

eraging over the different configurations (e.g. 2 Lc for a 
binary alloy) yields IV ,p(ii-u n ) and, by comparison with 
the free Green's function, the self-energy S(iuj n ). This 
fixes the l.h.s. of Eq. (|9"2"|) . i.e. the averaged Green's func- 
tion of the original model. Eq. ([92]) is then used to get 
a new IV ,p(i^ n ) and thus a new hybridization function 
via A(iw„) = iuj n + [i — t' — r) — G' {iu> n )~ l ■ New pa- 
rameters e[ X ^ and V^ X ^ arc obtained as the poles and 
weights of A.(iuj n ). From the above it is obvious, how- 
ever, that the self-consistency cycle can be set up for 
A directly. The variational determination of tg is not a 
problem: Analogous to the discussion given in Ref. loll 
a high-frequency expansion easily shows that, at station- 
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arity, t' must equal the intra-cluster part of the original 
hopping t. 

Note that for the non-interacting but disordered sys- 
tem, the M-CPA is related to the CPA as is, for the 
interacting but pure system, the cellular DMFT (C- 
DMFT)S£ to the DMFT. The M-CPA is a conceptually 
simple and straightforward generalization of the single- 
site CPA which includes short-range correlations in the 
disorder self-energy but is mean-field-like on a scale be- 
yond the cluster size. Like the variational cluster approx- 
imation (VCA), the M-CPA is systematic and controlled 
by 1/ L c as a small parameter, in principle. Practical cal- 
culations, however, are restricted to comparatively small 
clusters due to the exponential growth of the number 
of local (intra-cluster) configurations. For systems with 
local interactions and disorder, the structure of the ref- 
erence system Fig. [SJ: generates a combined C-DMFT 
+ M-CPA approach which is the straightforward clus- 
ter extension of the DMFT+CPA (Sec. [TXE)l . The ap- 
proximations within the SFT for different n s and L c are 
schematically grouped in Fig. [6] 



C. Translation symmetry 

There is an apparent problem remaining: Due to the 
construction of the reference system as a set of decou- 
pled clusters, the trial self-energies do not preserve the 
translational symmetries of the original lattice. This is a 
self-evident problem for any cluster approximation which 
is formulated in real space and has initiated the devel- 
opment of several further cluster approximations in the 
past. It has turned outfit however, that the construction 
of a self-consistent and systematic approach which on the 
one hand (for L c = 1) recovers the CPA and on the other 
(L c — ► oo) approaches the exact solution, which respects 
the requirements of causality and which perserves the 
lattice translational symmetries at the same time, is not 
easy to find. 

A straightforward idea is to use a reference system as 
displayed by Fig. 0). Here, isolated clusters are con- 
sidered again but with a hopping t' satisfying periodic 
boundary conditions. This restores translational invari- 
ance within the cluster at least. The idea has been put 
forward in the context of a pure but interacting system 
as a modified cluster-perturbation theory (CPT with pe- 
riodic boundary conditions) by Zacher et al4^ but was 
recognized^ to give less convincing results when com- 
pared to the usual CPT. Later on it could be showr*2£ 
within the VCA that periodic boundary conditions are 
in fact unfavorable: One simply has to treat the hopping 
integral connecting the edges of a cluster as a variational 
parameter. The SFT grand potential turns out to be 
stationary if this hopping integral vanishes. This corre- 
sponds to open boundary conditions. In any case, how- 
ever, it is obvious that the full translational symmetry 
cannot be restored with the reference system Fig. [TJd. 

Another straightforward idea is to formulate the orig- 



a) 




FIG. 7: (Color online) a) Anderson model, b) Possible ref- 
erence system consisting of isolated clusters with periodic 
boundary conditions generating a variant of the CPT. Within 
the VCA, the hopping marked as a dashed line may be treated 
as a variational parameter. 



inal problem as well as the reference system in recip- 
rocal space by using annihilators c& instead of c x etc. 
In graphical representations of H and H' like the one in 
Fig. dots would have to be reinterpreted as referring to 
one-particle states labeled by fc. The advantage is that, 
by construction, H' always exhibits the full translational 
symmetry. While there are no principle objections, it ap- 
pears to be impossible in practice, however, to generate 
meaningful approximations in this way. The reason is 
that a local interaction in real space transforms into a 
delocalized one in reciprocal space: The interaction pa- 
rameters Ukk'k"k" basically couple any k point to any 
other. Likewise, a local uncorrelated disorder transforms 
into a delocalized correlated one. As H and H' must 
share the same interaction and disorder, it is unlikely to 
find a t' that permits a (simple) solution of the resulting 
problem. 

A pragmatic way out would be to distinguish for- 
mally between the self-energy St>,p that is determined 
as the stationary point of the SFT grand potential on 
the one hand, and the translationally invariant ("phys- 
ical") self-energy on the other. The latter is obtained 
from S = Sf,p by some periodization procedure which 
employs a universal functional T "periodizing" the self- 
energy: S — ► T[S]. This idea has been suggested in the 
context of the C-DMFT42££ One possibility is to Fourier 
transform the optimized self-energy St\p from real space, 
S xx >, to reciprocal space via U Xi k, Eq. (|8T|) . This yields 
the self-energy in the representation Skk' ■ Substituting 
Skk' — * Skk^k.k' = T[S]kk' gives a translationally in- 
variant ("physical") self-energy T[S]. In real space, this 
periodization reads: 

TfS'jjcx' = ~j~ ^ ^ 3 X — X ' ,y — y f Syy' . (93) 

yy' 

T could also be applied to the Green's function calcu- 
lated from S via Dyson's equation. This is the usual 
procedure within the CPT to generate translationally in- 
variant Green's functions. While the problem of trans- 
lational symmetries can be fixed in this or in a similar 
way, the procedure appears to be ad hoc as it is placed on 
top of a variational (or self-consistent) calculation which 
itself involves S instead of T[S]. 
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D. Periodized M-CPA, periodized VCA 

Kollcr and Dupuis^i have discussed a modified form 
of the self-energy functional for (pure) systems of inter- 
acting bosons to get translationally invariant trial self- 
energies. This method can be adapted to disordered sys- 
tems as described below. 

Consider the Anderson model Eq. ([55)) and the follow- 
ing self-energy functional: 



n$,[s] 
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(94) 



As compared to the original functional f^pIS 1 ] given by 
Eq. (f64|) . the modified functional Sl^pfS'] differs in the 

first term as this includes the periodizing functional T. 
As usual, the hopping term of the Anderson model is sup- 
posed to be translationally invariant, T[t] = t. The exact 
self-energy is therefore translationally invariant, too, and 
satisfies T[St,p] = St.p- Hence, we have: 



'^p[s t , P ] = ii t ,p[St,p\ ■ 



(95) 



Using a one-particle basis of Bloch states labeled by wave 
vectors k to evaluate the traces, we furthermore have: 
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Sh',h - r[S]fe'fe . (96) 



Since T[St,p]kk' = 8k,k'F[St,p]kk, this shows that the ex- 
act (translationally invariant) self-energy is a stationary 
point of the modified self-energy functional. Conclud- 
ingly, the functionals fi£p[iS] and f^pfS] are different 
in general but coincide and are stationary at the exact 
self-energy St,p- They are similarly suited for construct- 
ing approximations. 

Consider now a trial self-energy St>,p from a reference 
system with the same disorder (i.e. the same distribu- 
tion function P) but with hopping parameters t' break- 
ing translational symmetries, as it is the case e.g. for the 
reference system given by Fig. [SJ;. The condition for sta- 
tionarity of the modified self-energy functional within the 
restricted set of trial self-energies then reads: 



= T 



n k,k' 
\T t ',p{iuJn) 



1 



iuj n + fi - t - T[S](iuj n ) r > 
dS kh ,{iw n ) 



k'k 



dt> 



(97) 



Here, S = St>,p for short. This replaces the condition 
Eq. (J90J) characteristic for the M-CPA. 

Fourier transformation U (Eq. (|87p) yields Skk' = 



E xx > e~ ikx e ik ' x ' S xxl and thus: 
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(98) 



For the reference system Fig. [5b with isolated clusters 
we have S xx > = S z+x . x '+x' = Ss+x,s+x'Sss' (and the 
same for its t! derivative) and thus the Euler equation 
(1981) is satisfied if 



ik(X-X') 



t iw„ + A*-e(fc)-T[S](fc) 
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(99) 



As compared to the self-consistency condition of the M- 
CPA, Eq. ([9^]) . the main difference consists in the fact 
that its solution is based on a self-consistency cycle which 
involves at each step the periodized self-energy T[S] in- 
stead of S. The Green's function of this periodized M- 
CPA, G = (iu n + i-i — t — TlS])^ 1 , is likewise translation- 
ally invariant: f[G]=G. 

The analog of the self-consistency equation (|99[) for the 
pure but interacting system is the self-consistency equa- 
tion of the periodized C-DMFT (PC-DMFT), see Ref. 
[50l . Hence, the above derivation also shows how the PC- 
DMFT can be rederived within the SFT framework. Fur- 
thermore, the construction of a combined PC-DMFT + 
periodized M-CPA is straightforward. It is also straight- 
forward to constuct a periodized VCA for disordered or 
for interacting systems along the lines above. 

As in the usual VCA (and also in the usual M-CPA) 
one in principle has the choice between open and peri- 
odic boundary conditions in the reference system. In the 
case of the periodized VCA (periodized M-CPA) a trans- 
lationally invariant self-energy and Green's function are 
generated by either choice. The value of the SFT grand 
potential at the stationary point decides which kind of 
boundary conditions should be preferred. 

With a suitably defined functional T there are no prob- 
lems to extend the approach and to restore point-group 
symmetries of the underlying lattice which could be vi- 
olated by the choice of the reference system. Essentially 
the same arguments given for the case of translational 
symmetries can be repeated. 

There is yet another way to modify the original func- 
tional and thereby to construct a theory which respects 
the translational symmetries of the underlying lattice. 
This shall be mentioned here for the sake of complete- 
ness. Consider the functional 
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Again, only the first term on the r.h.s. is modified. Note 
that this is essential for the construction of approxima- 
tions: It ensures that the not explicitly known but uni- 
versal functional given by the remaining terms cancels 
out when comparing with the functional of a suitably de- 
fined reference system, as usual, P [S] coincides with 

f2t,p[S] and Sl^pfS'] at the exact self-energy St,p and is 
stationary there. With the reference system Fig. [5t the 
following self-consistency equation can be derived: 



}_^ e ik(X-X') 
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(101) 

As compared to the periodized M-CPA, the conceptual 
disadvantage consists in the fact that only the propagator 
T [(GtQ — S') -1 ] but not the self-energy is translation- 
ally invariant. 



E. Dynamical cluster approximation 

Originally, the dynamical cluster approximation 
(DCA) has been proposed as a cluster extension of the 
DMFT for interacting systems4^££ Essentially the 
same ideas, however, can also be used to construct a 
generalization of the single-site CPA for the disorder 
problemi 44 i 45 Here, it is shown that the (disorder) DCA 
can be rederived within the SFT framework by utilizing 
the real-space perspective on the DCA first discussed by 
Biroli et al^ 

While the main idea of the periodized M-CPA to re- 
store translational symmetry is to consider a modified 
but equivalent self-energy functional, one could also keep 
the exact functional form f^^S 1 ] but modify the hop- 
ping of the original system, i.e. t — ► t. Approximations 
are then constructed by starting from ilj P [S] and using a 
reference system consisting of isolated clusters again. To 
ensure that the resulting approximations systematically 
approach the exact solution for cluster size L c — > oo, the 
replacement t — ► t must be controlled by L c , i.e. it must 
become exact (up to irrelavant boundary terms) in the 
infinite-cluster limit. 

Consider, for example, 



t = (VW)UHU(VW) 



(102) 



For clusters of finite size L c , the combined Fourier trans- 
formation VW is different from U. For L c — > oo, how- 
ever, this becomes irrelevant. With e(k) = (UHlJ)(k) 
we have: 
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e(k + K) 



(103) 

Obviously, t is invariant under superlattice translations 
as well as under cluster translations (with periodic clus- 
ter boundary conditions). The original and the modified 
system with H = H(t,r],0) are represented by Fig. [5^,, 



b. The construction of t is such that it exhibits the same 
translational symmetries as the one-particle parameters 
t' of a reference system consisting of isolated clusters 
tiling the original lattice with periodic boundary condi- 
tions, see Fig. [8p, d. Since both, t and t, are invariant un- 
der superlattice translations, we can compare txx 1 (fe) = 

(vHv) xx >(k) with t xx ,{%) = (vHv) xx ,(k). it 

turns out they are equal up to a phase factor: 
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(104) 



To rederive the disorder DCA within the framework of 
the SFT, we use the functional flj p [S] and the reference 
system Fig. [SJi. Therewith, one formally arrives at the 
M-CPA self-consistency condition, Eq. (|9"2"|) . but with t 
replaced by t: 
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(105) 

The bold symbols are matrices in the cluster variables 
X,X'. The self-energy and the averaged Green's func- 
tion of the reference system are independent of k. The 
decisive difference as compared to the M-CPA is that the 
(modified) hopping of the original system t is invariant 
under cluster translations, as it is the case for t' . This is 
important as it allows to simultaneously diagonalize all 
matrices in Eq. (|105j) by the cluster Fourier transforma- 
tion W . Furthermore, from the definition of t, the W 
transformation yields t(k) — > e(k + K) = e(k), i.e. we 
find: 
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This, however, is just the self-consistency equation of the 
disorder DCA4i 

We can thus state that, analogous to the work of Biroli 
et ali££ for the usual DCA, it is found that the disor- 
der DCA is equivalent with the M-CPA applied to the 
system with modified hopping Eq. (|102p . Hence, it has 
been shown that the (disorder) DCA (as well as the usual 
DCA and also the combined theory for interacting sys- 
tems with disorder) can be recovered within the frame- 
work of the SFT when starting from a suitably modified 
problem. As already been noticed,— however, a strict 
rederivation starting from the original system (with t) 
appears to be impossible. 

In this context an interesting new approximation sug- 
gests itself: Starting with the modified hopping t and 
using the reference system shown in Fig. [5Ji, generates a 
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FIG. 8: (Color online) a) Anderson model, b) The original 
Anderson model but with a modified one-particle part t — ► t 
which is the starting point for the dynamical cluster approx- 
imation (DC A), t is invariant under superlattice and cluster 
translations, c) Reference system generating the DCA. Note 
that t' has the same translational symmetries as t. d) Refer- 
ence system generating a simplified DCA (see text). 



simplified (disorder) DCA without bath degrees of free- 
dom. This simplfied DCA is a systematic (controlled 
by l/L c ) cluster approximation and gives a translation- 
ally invariant self-energy and propagator. The simplified 
DCA is related to the periodized VCA in the same way 
as the DCA is related to the periodized M-CPA. The 
analogous formulation of a simplified DCA for the pure 
but interacting system represents a variational extension 
of a non-self-consistent approximation ( "periodic CPT" ) 
recently introduced by Minh-Tien4£ 



XI. SUMMARY AND DISCUSSION 

An advantageous feature of the self-energy-functional 
approach for systems of interacting electrons is that ap- 
proximations are easily specified by choosing a reference 
system. The reference system helps to span a space of 
trial self-energies which are optimized using an appro- 
priate dynamical variational principle. Since a reference 
system must share with the original system under con- 
sideration the same interaction part, the number of pos- 
sible reference systems is severly limited. This leads to 
a straightforward classification of approximations which 
may be called "dynamic" as these are essentially based 
on an approximation for the self-energy or, equivalently, 
the Green's function, i.e. quantities characterizing the 
spectrum of one-particle excitations. 

For systems with local interactions, the set of dynam- 
ical approximations includes the dynamical mean-field 
theory which corresponds to a reference system consist- 
ing of decoupled single-impurity models. Modifying this 
reference system generates DMFT-related dynamical ap- 
proximations but does not spoil the main attractive prop- 
erties of the DMFT, in particular its non-perturbativc 
character and its thermodynamical consistency. This is 
interesting as it opens up a way to construct approxima- 
tions which (i) are based on simpler reference systems 
that can be solved with less numerical effort or/and (ii) 
include spatial correlations beyond the local mean-field 



concept. In fact, the classification of dynamical approxi- 
mations has led to new approximations which have suc- 
cessfully been employed in the past (see Fig. [5] for pure 
systems) : The dynamical- impurity approximation (DIA) 
as well as the variational cluster approximation (VCA). 

The main idea of this paper has been to translate this 
story to the case of disordered (and interacting) systems. 
The disorder is assumed to be local and uncorrelated be- 
tween different sites. One possible way is to apply the 
conventional SFT to treat the interaction part of the 
problem for any disorder configuration and to average 
subsequently. This procedure treats the disorder part 
of the problem exactly. Choosing a decoupled set of ef- 
fective impurity models as a reference system, yields an 
extension of the DMFT for systems with reduced transla- 
tional symmetries or, after reinterpretation of the mean- 
field (Euler) equations as stochastic recursion equations, 
to the so-called statistical DMFT. This is feasible in the 
case of a Bethe lattice only. Even then and even with 
a simpler reference system including a minimum num- 
ber of bath sites, however, the statistical SFT remains a 
numerically extremely expensive method. 

The alternative consists in an approximate treatment 
of the disorder on the same footing as the interaction. 
Technically, this requires a reformulation of the SFT 
for interacting systems within a functional-integral lan- 
guage which, as a by-product, provides an entirely non- 
perturbative construction of the SFT, i.e. avoids for- 
mal summations of skeleton diagrams. The functional- 
integral framework then allows to formulate a disorder 
SFT in essentially the same way as for pure interact- 
ing systems; the main corners of the theory are left 
unchanged: The (averaged) grand potential can be ex- 
pressed as a functional of the (disorder) self-energy with 
the physical self-energy being a stationary point. The 
functional still consists of a simple and explicitly known 
part depending on the one-particle parameters as well as 
of a complicated and basically unknown part which, how- 
ever, is "universal" , i.e. depends on the probability distri- 
bution only. Choosing a reference system with the same 
distribution function, the universal part can be elimi- 
nated, and an exact evaluation of the self-energy func- 
tional becomes possible on the space of trial self-energies 
generated by the reference system. Restricting the search 
for the stationary point to this limited subspace gener- 
ates approximations. 

It is conceptually very satisfying that within this gen- 
eralized SFT, the coherent-potential approximation for 
the disorder problem takes the place of the DMFT for 
an interacting system. As the DMFT, the CPA is a dy- 
namical approximation for the self-energy and is distin- 
guished by the fact that it becomes formally exact in the 
limit of infinite spatial dimensions. "Replacing interac- 
tion by disorder" , any reference system for the case of an 
interacting system can also be considered as a reference 
system for the disorder case. This establishes a one-to- 
one mapping of the respective approximations with the 
DMFT corresponding to the CPA and the cellular DMFT 
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corresponding to the molecular CPA. Via this mapping 
a disorder DIA and a disorder VCA appear as new ap- 
proximations as well as the periodized M-CPA, namely 
the disorder pendant of the periodized cellular DMFT, 
when starting from a different but equivalent functional 
involving the periodizing functional T. Starting from 
the self-energy functional with a suitably modified hop- 
ping t — ► i, one also recovers the disorder analog of the 
dynamical cluster approximation. Finally, a simplified 
DCA without bath degrees of freedom can be set up. 

The case of non-interacting disordered systems could 
actually be treated by specialization of a more general 
functional Qt,p,u[S, {^ij}] depending on the full and 
on the configuration-dependent interaction self-energies. 
This is applicable to interacting and disorderd systems. 
Though formally more complicated, the not explicitly 
known part of the functional is universal, i.e. depends 
on U and P only. The classification of dynamical ap- 
proximations extends accordingly (sec Fig.[6|). 

Summing up, the type of dynamical approximations 
that can be constructed for systems with local interac- 
tions or/and local disorder are mean-field and cluster ap- 
proximations which differ with respect to the number of 
local variational degrees of freedom included in the refer- 
ence system. All approximations fulfill the requirements 
of causality. They are non-pcrturbative as, contrary to 
truncating diagrammtic approximations, the exact func- 
tional form is retained. Thermodynamical consistency is 
ensured by the existence of an approximate but explicit 
expression for the grand potential, i.e. for a thermody- 
namical potential, from which the physical quantities are 
derived. Additional but systematic modifications of the 
sclf-encrgy functional have been shown to generate (clus- 
ter) approximations that respect the symmetries of the 
underlying lattice, particularly the translational symme- 
tries. 

Convergence properties of the different quantum- 
cluster schemes with a continuum of bath sites have been 
discussed in Refs. [88. 89. 90. 9l| and apply to the schemes 
with n s = 1 in an analogous way. For L c — > oo, local 
quantities generally converge exponentially fast within 
the C-DMFT / M-CPA and within the (disorder) VCA. 
This has to be compared with the l/l^ behavior obtained 
within the DCA and within the simplified DCA (with 
L c = I®). For practical purposes, however, the quality 
of a given approximation for small L c is more important 
and can apparently be estimated a posteriori only. It is, 
for example, an open question whether cluster schemes 
with or without bath degrees of freedom should be pre- 
ferred in this respect (see, however, Ref. 1281). 

The application of the SFT has been restricted to sys- 
tems with local interactions and local disorder. This is 
consistent with the spirit of the cluster-mean-field ap- 
proximations discussed above. One should note, however, 
that non-local or even long-ranged interactions or disor- 
dered nearest-neighbor or longer-ranged hopping or spa- 
tially correlated disorder in the on-site energies, for exam- 
ple, pose difficulties. While the general self-energy func- 



tional can be set up as usual, it appears to be impossible 
to find a suitable reference system and thereby usable 
approximations in most cases: As the reference system 
should have the same interaction and disorder, it is hard 
to find a decoupling of its degrees of freedom by modify- 
ing the hopping part only. This problem could be han- 
dled either pragmatically by additional mean-field decou- 
plings of non-local terms connecting different clusters^ 
or, more thoroughly by considering more complex func- 
tionals involving two-particle correlation functions^ 

One should be aware that the various local approx- 
imations and their different cluster extensions all ne- 
glect long-range correlations beyond the linear scale of 
the cluster. This is typical for any cluster mean-field 
approach. To our knowledge there is no possibility to 
systematically restore long-range correlations by a suit- 
able embedding of an isolated cluster in the environment 
within the presented formalism (apart from simply en- 
larging the cluster). The same holds for cluster schemes 
formulated within reciprocal space. 

This implies that a proper description of transport 
properties is hardly possible. Effects like Anderson local- 
ization, for example, cannot be captured within a mean- 
field approach like CPA and are obviously difficult to 
restore by cluster schemes extending CPA (see Ref. H . 
for an example). Anyway, it is basically impossible to 
access two-particle correlation functions, and the con- 
ductivity in particular, within an approach that places 
one-partible excitations in the center of interest. The sta- 
tistical DMFT (statistical SFT) represents an exception 
as here the full distribution of the (one-particle) density 
of states can be used to discriminate between extended 
and localized states, for example. 

Besides the spectrum of one-particle excitations, how- 
ever, the SFT derives an approximate but explicit ex- 
pression for a thermodynamical potential and thus pro- 
vides a consistent picture of the entire thermodynamics 
and of static expectation values. This includes sponta- 
neous symmetry breaking, i.e. the determination of or- 
der parameters. The formal framework presented here 
should therefore be ideally suited to study the effects 
of disorder on different types of long-range order, such 
as magnetism or superconductivity. For different ma- 
terial classes, such as diluted magnetic semiconductors, 
cuprate-based high-temperature superconductors, man- 
ganites, rare-earth compounds etc. these are central ques- 
tions. 
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APPENDIX A: DISORDER IN THE 
INTERACTION 

It is formally straightforward to extend the theory to 
the case of disorder in the interaction part of the Hamil- 
tonian. The theory can be constructed without major 
modifications: 

Consider a system with Hamiltonian H = Ho(t) + 
Hdis(v) + H'mt(U) where H (t) and H- lnt {U) describe the 
one-particle and the interaction part as before while 



t t 



(Al) 



is an interaction term with random parameters r/ dis- 
tributed according to some P{ri). The definition of the 
free Green's function Eq. (fTT|) . the interacting Green's 
function Eq. (fl~3|) . and the averaged Green's function Eq. 
(|16p is formally unchanged. The self-energies are defined 
as 

^t, v ,u = G t o o - G t l l U , S t ,p.u = G t o o - r t PE/ . 

(A2) 

The reasoning in section Hill is unchanged. The function- 
als [G ■ ^u.tiIGq 1 ], Gu.n [G , Gf/^P], 

and Fu^[S\, however, additionally depend on the pa- 
rameters rj. The same holds for section I VII except for 
Eq. W!\ which has to be replaced by 
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The final self-energy functional thus reads 
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Approximations are constructed in the same way as de- 
scribed in Sec. I VIII by making contact with a reference 
system with the same interaction and disorder. 



APPENDIX B: THERMODYNAMIC AL 
CONSISTENCY 



particle number 

Nt,p,u = J drjP{r])N w = J dr)P{ri)tx{ PtiVtU N) 

(Bl) 

can be obtained: (i) Nt,p : u is calculated, on the "zero- 
particle level" as the fi derivative of the averaged grand 
potential: 

dflt,p,u 



t.P.U 
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(B2) 



or (ii) N t ,p,u is calculated, on the "one-particle level" 
by frequency integration of the one-particle averaged re- 
tarded Green's function: 

Nt.P.u = / duf(uj)Jmtrr tt p tU (uj + t0+) , (B3) 

7T J-oo 

where f{w) = l/(e"/ T + 1) is the Fermi function. Ther- 
modynamical consistency means that both ways must 
yield the same result. This is not clear a priori as within 
the SFT Ctt.p.u and ^t,p,u ar e approximate quantities. 

As in Ref. [61], however, it can be argued that there is 
a twofold fi dependence of Qt,p.u- an explicit /i depen- 
dence which in Eq. ((50]) is due to the free Green's function 
Gto o = u + l L — *j an( i an implicit fj, dependence due to 
the fjj dependence of the disorder self-energy at the sta- 
tionary point. Now, if and only if an overall shift s of 
the on-site energies in the reference system is treated as 
a variational parameter, the derivative w.r.t. the implicit 
(i dependence vanishes, and one is left with the explicit 
one: 



t,p,u 



t,P.U 



da 



t,p,u 



dfx 




dfi, 



cxpl. 



C^oxpl. 



Trln- 



Moxpl 



- 1 - S t > 



v,p,u 



(B4) 



according to Eq. (|60|) . Carrying out the differentiation 
and using ^t,p,u = (G^q — St' .p.u)^ 1 , one immediately 
arrives at Eq. (|B3[) which proves the equivalence with Eq. 
(|B2|). 

Typically, only a distinguished set of parameters are 
treated as variational parameters in a practical calcula- 
tion. The argument shows the necessity to include e' 
in the set of variational parameters if thermodynamical 
consistency is required for the particle number. The ar- 
gument straightforwardly generalizes to all one-particle 
operators coupling linearly with a parameter A' to the 
reference-system Hamiltonian. 



Here, we briefly discuss the thermodynamical consis- 
tency of approximations within the SFT generalized to 
disordered systems. To be definite, we concentrate on the 
particle number as an example. The reasoning closely 
follows Ref. H 

There are two ways by which the configurational av- 
erage of the quantum-statistical expectation value of the 



APPENDIX C: CAUSALITY 

The causality of all (approximate) dynamical quanti- 
ties (configuration-dependent as well as averaged Green's 
functions and configuration-dependent/independent self- 
energies) is inevitable for having well-defined self-energy 
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functionals and expressions involving Trln(---) in par- 
ticular. A frequency-dependent quantity is termed to be 
causal if it is analytical in the entire complex uj plane 
except for first-order poles on the real axis with posi- 
tive residues. Equivalently, we can demand that, after 
retardation uj — > uj + i0 + (u) real), the imaginary part be 
negative definite. 

Causality is easily verified, for example, for all quanti- 
ties which appear in the self-energy functional Eq. (|5Tj) 
when evaluated on its domain (see Sec. IVIIIjl . The last 
three terms on the r.h.s. only involve exact and therewith 
causal quantities as arguments of Trln(- • • ) (as these are 
taken from the reference system). Note that the configu- 
ration average of a configuration-dependent exact quan- 
tity is causal: 

Let G(uj) be the exact Green's function of a model 
H. We have the Lehmann rcprcsentationiS G(uj) = 
Qg{uj)Q^ with ^-independent matrices Q (Q not 
quadratic, see Ref. l93l ) and a diagonal matrix g(uj) 
with elements g nn (uj) = l/C^ — Wn), where uj n are 
the poles of G. Consider T(uj) = (G(uj)} averaged 
over only two configurations, for simplicity, i.e. consider 
T(ui) = PiQigi(u)Q\ + P2Q292MQ2 with pi,p 2 > 



0,Pi +P2 = 1- With Q = (y/piQi, and j(u>) 

[ Q 1 ^ ' ? , I we then immediately have the repre- 
V 92 (w) J J 

sentation T(u>) = Qj(uj)Q^ from which the causality of 

r(w) is easily verified. 

The only possible source of causality violation is left 
for the first term on the r.h.s. of Eq. (|5ip . However, it 
is easily shown that (i) a causal T'(ui) (the exact aver- 
aged Green's function of the reference system) implies 
a causal S(uj). The main point is that the retarded 
quantity S(ui + i0 + ) = uj + n - t' - T'(uj + iO+y 1 has 
a negative definite imaginary part since the imaginary 
part of T'(lu + iQ + )~ 1 is positive definite according to 
a lemma given in Ref. I27I . Furthermore (ii), it can be 
shown that the causality of S(uj) implies the causal- 
ity of r(w) = (Go(w)" 1 - S(w))" 1 for arbitrary G = 
[io + ^-t)- 1 . Namely, with S{ui) = G' {uj)- 1 -V {uj)- 1 
and V(u>) = Q'^'{u>)Q we have the representation 
T{ui) = l/((Q'7 , (w)Q' t )- 1 - (t - t')) and thus (see Ref. 

M) r(w) = Q / (7'(w)- 1 -Q /t (i-t')QT 1 Q /t from which 
is it obvious that T(uj) is analytical with the exception 
of first-order real poles with positive residues. 
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